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ABSTRACT 


Target  recognition  and  identification  in  battlefields  has  been  a  crucial  determinant 
to  the  ultimate  success  or  failure  of  modern  military  campaigns.  Since  World  War  II,  the 
Identification  of  Friend  or  Foe  (IFF)  systems  installed  in  radar  systems  have  served  as  the 
primary  cooperative  target  identification  techniques  based  on  the  “question  and  answer” 
interrogation  loop  of  unidentified  aircraft.  However,  the  IFF  system  has  a  number  of 
limitations  that  pose  significant  challenges  in  the  positive  identification  of  hostile  and 
neutral  aircrafts,  which  can  lead  to  a  catastrophic  outcome  of  fratricide  and  the  possible 
elimination  of  a  friendly  or  commercial  aircraft. 

To  reduce  the  probability  of  fratricide  and  enhance  the  effectiveness  and 
robustness  of  target  identification,  this  research  will  examine  the  radar  Non-Cooperative 
Target  Recognition  (NCTR)  technique  of  using  the  bi-spectral  signatures  of  backscattered 
radar  signals.  The  basic  idea  is  that  the  geometry  of  the  target  scatterers  and  their  mutual 
interactions  impose  features  in  the  reflected  radar  signal  that  are  typical  and  unique  to  the 
target  of  interests.  The  bi-spectrum  can  be  used  to  detect  these  multiple  interactions 
features  which  may  then  be  used  to  match  against  a  reference  database  that  contains 
signatures  of  different  target  types  for  recognition  and  identification. 
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I.  INTRODUCTION 


Know  the  enemy  and  know  yourself;  in  a  hundred  battles  you  will  never 
be  in  peril.  Know  yourself  but  not  the  enemy;  for  every  victory  gained 
you  will  also  suffer  a  defeat.  Know  neither  the  enemy  nor  yourself;  you 
will  succumb  in  every  battle. 

—  Sun  Tzu,  Chinese  General  and  Military  Strategist  (1971) 


A.  MOTIVATION 

Situation  awareness  has  always  been  the  key  determinant  for  the  ultimate  victory 
or  defeat  of  all  warfare  and  military  campaigns.  Through  situation  awareness,  one  must 
be  able  to  ascertain  the  four  “Ws”  of  one’s  own  forces  as  well  as  one’s  enemies,  which 
are: 

1 .  Where  are  they:  azimuth,  elevation,  range,  etc? 

2.  When  will  they  arrive:  heading,  velocity,  time  of  flight,  etc? 

3.  Who  are  they:  friendly,  hostile  or  neutral?  and 

4.  What  are  they:  fighter  jet,  commercial  aircraft,  tank,  school  bus,  etc? 

While  “Where”  and  “When”  are  adequately  addressed  by  the  early  RAdio 
Detection  And  Ranging  (Radar)  systems  used  primarily  as  sensors  to  detect  and  track 
targets  since  World  War  II,  “Who”  and  “What”  or  more  specifically,  target  recognition 
and  identification,  can  pose  significant  challenges,  as  only  limited  information  pertaining 
to  the  characteristics  of  the  target  being  detected  has  been  available  from  the  radar 
backscattered  signal. 

As  such,  early  radar  systems  were  merely  “blob”  detectors,  as  detections  of  target 
presences  from  the  reflected  radar  signal  were  usually  displayed  as  phosphor  “blips”  on 
the  radar  system’s  Plan  Position  Indicator  (PPI)  (Skolnik,  2001).  Information  on  the 
target’s  location  in  range,  azimuth  and  elevation  and  velocity  could  then  be  extracted 
from  the  PPI,  but  could  not  provide  further  detail  on  the  target’s  type  and  identity. 
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This  limitation  of  the  early  radar  systems,  together  with  the  vital  need  for  situation 
awareness  as  a  force  multiplier  to  enhance  combat  effectiveness,  provided  the  impetus  for 
the  development  of  the  IFF  system  and  various  other  radar  target  recognition  techniques. 

1.  IFF  System  and  Cooperative  Target  Recognition 

The  IFF  system  is  a  two-channel  interrogator  /  transponder  system,  with  one 
frequency  used  for  the  interrogating  signals  and  another  for  the  response,  that  uses  pulse- 
coded  radio  frequency  transmissions  to  positively  identify  friendly  aircraft,  ships,  and 
army  units  from  those  of  the  enemy  forces.  It  serves  as  a  cooperative  technique  for  target 
recognition,  based  on  the  “Question”  interrogation  and  “Answer”  response  loop  of 
unidentified  targets  of  interests,  as  depicted  in  Figure  1 . 


Figure  1  IFF  “Question  and  Answer”  Loop  (Wolff,  1997) 


The  IFF  interrogator  system  initiates  the  “Question  and  Answer”  loop  by 
transmitting  a  pulse-coded  radio  frequency  signal,  also  referred  to  as  the  “challenge”  that 
asks  a  series  of  questions,  one  of  which  is  “Who  are  you?”  to  the  unidentified  target  of 
interest.  The  interrogation  code  or  challenge  is  then  received  by  an  electronic  system 
known  as  a  transponder  that  is  installed  onboard  the  target.  If  the  target  is  friendly  and 
equipped  with  a  compatible  transponder,  it  automatically  interprets  the  challenge  and 
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responds  by  transmitting  the  requested  identification  code  back  to  the  interrogator 
system.  After  successfully  receiving  and  decoding  the  acknowledgement,  the  target  is 
tagged  as  a  friendly  target  by  the  interrogating  radar  system. 

Conversely,  if  the  target  is  neutral  or  hostile,  it  is  assumed  that  the  challenge 
cannot  be  decoded  and  interpreted  by  the  onboard  transponder.  Henceforth,  the 
transponder  is  unable  to  respond  with  a  correct  answer  and  will  be  tagged  as  an  unknown 
or  foe  target  by  the  interrogating  radar  system. 

Although  the  IFF  system  plays  a  significant  role  in  contributing  toward  target 
recognition,  it  can  only  perform  positive  recognition  of  friendly  targets,  and  it  requires 
the  active  and  cooperative  participation  of  the  targets  which  are  under  interrogation. 

Moreover,  IFF  systems  have  a  number  of  limitations,  and  cooperative  systems 
have  historically  been  plagued  with  the  following  fundamental  issues: 

a.  If  the  IFF  system  is  inoperative,  due  to  a  malfunction,  or  an  operator 
intentionally  turns  it  off,  there  is  no  inherent  method  for  a  weapon  system  to 
distinguish  the  resulting  lack  of  response  from  that  of  an  enemy  target. 

b.  Hostile  targets  may  intercept  /  listen  to  and  re-transmit  an 
acknowledgement  from  a  “friendly”  target,  in  which  case  the  hostile  target  will  be 
tagged  as  “friendly.” 

c.  Commercial  airliners  are  fitted  with  transponders  that  conform  to  open 
international  civil  aviation  standards,  which  enable  terrorists  or  enemy  forces  the 
opportunity  to  be  equipped  with  such  transponders  on  their  aircrafts. 
Consequently,  they  can  masquerade  as  passenger  airliners  for  the  purpose  of 
infiltrating  the  air  defence  system. 

Hence,  it  is  obvious  that  the  IFF  system  is  not  infallible  and  any  interruption  in 
the  IFF  “Question  and  Answer”  loop,  resulting  in  the  absence  of  a  response  or  correct 
acknowledgement,  can  potentially  draw  the  inference  that  the  unidentified  target  is 
hostile.  This  seriously  reduces  the  utility  of  neutral  or  allied  forces  and  can  potentially 
lead  to  catastrophic  outcomes  of  fratricide  and  the  erroneous  elimination  of  a  friendly  or 
commercial  aircraft. 
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2.  Tragic  Incidents  of  Fratricide 

To  illustrate  the  devastating  consequences  of  erroneous  target  recognition  and 
identification,  a  few  examples  of  wrongful  engagements  of  neutral  civilian  and  friendly 
targets  from  recent  history  are  referenced  from  V.G  Nebabin’s  Methods  and  Techniques 
of  Radar  Recognition  (1995)  and  P.  Tab’s  Introduction  to  Radar  Target  Recognition 
(2005). 

a.  In  1980,  a  DC-9  Italian  passenger  plane  was  shot  down  by  mistake  with  a 
missile  (supposedly  American  or  French  made)  near  the  Island  of  Sicily.  It  is 
believed  that  the  missile  was  re-aimed  by  mistake  because  of  incorrect  IFF  and 
discrimination  of  the  type  of  plane.  Eighty-one  passengers  died. 

b.  In  1983,  a  Boeing  747  passenger  plane,  South  Korean  Airlines  Flight  007, 
flying  from  New  York  and  Anchorage  to  Korea,  was  shot  down  by  the  USSR  Air 
Defense  SU-15  fighter  as  a  result  of  incorrect  recognition.  The  Boeing  747  was 
recognized  as  an  RC-135  reconnaissance  aircraft.  Two-hundred  and  sixty-nine 
passengers  from  Korea,  Japan,  the  U.S.  and  other  countries  died. 

c.  In  1988,  an  Iranian  A-300B  Airbus  passenger  plane,  Flight  655,  flying 
from  Bender-Abba,  Iran,  to  Dubai,  United  Arab  Emirates  was  shot  down  by  an 
SM-2  Standard  Missile  under  control  of  the  Aegis  system  of  the  cruiser  USS 
Vincennes  in  the  Persian  Gulf.  The  A-300B  was  recognized  erroneously  as  an  F- 
14  fighter.  Two-hundred  and  ninety-eight  people  died. 

d.  In  1994,  two  U.S.  Blackhawk  military  helicopters  were  shot  down  in  the 
north  of  Iraq  by  two  U.S.  F-15  fighters.  The  passengers  and  crews  of  the 
helicopters  had  been  performing  a  humanitarian  mission.  They  had  been 
recognized  erroneously  as  Iraqi  helicopters  flying  in  the  no-fly  zone.  Fourteen 
men  died. 

e.  During  the  2003  Iraq  war,  there  were  several  incidents  of  ground  troops 
being  killed  by  friendly  fire  from  the  air.  A  UK  Tornado  aircraft  was  shot  down 
by  a  Patriot  missile,  apparently  being  mistaken  for  a  Scud  ballistic  missile,  and 
two  British  soldiers  were  killed  apparently  in  an  exchange  of  fire  between  two 
British  Challenger  tanks. 
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3.  Non-Cooperative  Target  Recognition  (NCTR) 

To  prevent  the  catastrophic  outcome  of  fratricide  and  enhance  the  robustness  of 
the  air  defense  system,  NCTR,  which  encompasses  the  techniques  and  technologies  to 
extract  and  examine  the  characteristics  of  target  signatures  in  order  to  provide  further 
insight  on  the  target’s  type  and  identity  without  its  active  cooperative  participation,  has 
been  developed  and  adopted  over  the  years  (Tait,  2005). 

A  variety  of  active  and  passive  sensor  systems  can  be  utilized  to  perform  NCTR 
and  establish  a  target’s  type  and  identity.  A  discussion  of  the  various  sensor  systems  and 
their  strengths  and  limitations  in  target  recognition  perfonnance  is  presented  in  the 
following  sections. 

a.  Electro-Optical  (EO)  and  Infrared  (IR)  Systems 

An  EO  system  uses  optics  for  the  surveillance,  tracking  and  recognition  of 
a  target  and  depends  on  external  illumination  for  the  conversion  of  photons  into  amplifier 
electrical  signals  for  imaging.  An  IR  system  detects  the  heat  emission  of  the  target  and 
can  operate  in  total  darkness  without  external  illumination  of  the  target  as  it  is  dependent 
on  the  temperature  differences  between  the  target  and  its  background,  for  imaging.  Both 
EO  and  IR  systems  are  passive  sensor  systems  that  are  ideal  for  covert  operations  and 
have  strength  in  providing  high-resolution  images  for  target  recognition.  However,  their 
performances  are  adversely  affected  by  weather  effects,  such  as  rain,  fog  and  haze,  and 
by  atmospheric  attenuation  due  to  scattering  and  absorption. 

b.  Electronic  Warfare  (EW)  and  Electronic  Support  Measure 
(ESM)  Systems 

EW  and  ESM  systems  use  passive  receivers  for  the  detection  and  analysis 
of  signal  in  the  electromagnetic  spectrum.  The  characteristics  of  the  intercept  signal  are 
measured  and  used  to  isolate,  identify  and  locate  the  emitter.  Identification  of  the  emitter 
through  the  use  of  finger-printing  techniques  of  its  measured  signal,  along  with  a  library 
of  reference  signals  can  then  be  used  to  recognize  the  associate  platform.  Although  EW 
and  ESM  systems  have  good  range  perfonnance,  due  to  the  one-way  propagation  of  the 
emitter  signal,  their  effectiveness  is  completely  reliant  on  the  target’s  electromagnetic 
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emission  for  locating  and  possibly  identifying  it.  Hence  an  Emission  Control  (EMCON) 
imposed  by  an  incoming  target  will  easily  overcome  the  EW  and  ESM  systems’  NCTR 
capability. 

c.  Acoustics  System 

The  acoustics  system  uses  an  array  of  microphones  to  discern  sound 
waves  from  the  noisy  environment  for  the  detection,  localization  and  possible  recognition 
of  sound  generating  targets,  such  as  human  speech,  the  engine  of  a  moving  vehicle  and 
the  chopping  blades  of  a  hovering  helicopter,  etc.  While  the  acoustics  signal  from  these 
targets  may  be  discriminating  and  recognisable,  the  acoustics  sensor  is  generally  limited 
in  performance  due  to  the  slow  speed  of  sound,  which  is  about  330  m/s  in  air,  and  limited 
sensitivity  to  background  noise.  The  acoustics  system  is  also  affected  by  atmospheric 
propagation  attenuation  and  is  only  considered  for  very  short-range  applications. 

d.  Radar  System 

By  far,  a  radar  system  is  considered  the  sensor  of  choice  for  NCTR,  as  it 
offers  the  most  potential.  It  is  able  to  perform  round-the-clock  24/7  surveillance,  can 
perform  detection  and  tracking  of  targets  at  ranges  much  further  than  EO  and  IR  systems, 
and  is  not  affected  by  adverse  weather  conditions.  Furthermore,  unlike  EW  and  ESW 
systems,  radar  systems  are  not  dependent  on  the  target  emissions.  Technological 
advancements  have  created  a  system  that  formerly  could  only  detect  and  track  targets  to 
one  that  is  capable  of  high  resolution  imaging  and  target  recognition. 

According  to  Van  Der  Heigen’s  “Aircraft  Recognition  with  Radar  Range 
Profiles  (1998),”  Radar  NCTR  techniques  are  rooted  on  the  basic  concept  that  the 
geometry  of  the  target  or  the  moving  parts  on  the  target  impose  characteristic  features  in 
the  reflected  radar  signal  that  are  typical  and  are  unique  to  the  targets  of  interest.  These 
features  may  then  be  used  to  match  against  a  reference  database  that  contains  signatures 
of  different  target  types  for  target  recognition  and  identification. 
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B.  SCOPE  AND  ORGANIZATION 

1.  Thesis  Scope 

This  thesis  first  discusses  the  two  fundamental  methods  of  High  Range  Resolution 
(HRR)  imaging  and  Signal  Modulation  analysis,  developed  to  perform  target  recognition 
and  identification  using  a  radar  system.  Several  conventional  NCTR  techniques  based  on 
these  two  fundamental  methods  are  then  reviewed  and  evaluated  based  on  their  relative 
merits  and  limitations. 

Next,  the  thesis  focuses  on  developing  a  complementary  radar  NCTR  technique  of 
using  the  bi-spectral  signatures  of  backscattered  radar  signals  to  augment  the  knowledge 
that  conventional  NCTR  techniques  already  offer.  Derivation  and  theoretical  calculations 
of  the  bi-spectrum  of  radar  target  signature  were  performed  to  extract  the  multiple 
scattering  interaction  features  which  may  then  be  used  to  match  against  a  reference 
database  that  contains  signatures  of  different  target  types  for  target  recognition  and 
identification. 

2.  Research  Objectives 

The  objectives  of  this  thesis  are  two-fold: 

a.  To  study  and  mathematically  model  the  direct  and  indirect  multiple 
interactions  scattering  mechanisms  of  radar  backscattered  signal. 

b.  To  develop  the  radar  NCTR  technique  of  using  the  bi-spectral 
signatures  of  backscattered  radar  signals  for  characterization  of  a  target’s  multiple 
interactions  features. 

3.  Outline  and  Organization 

Chapter  II  provides  a  review  of  the  current  conventional  NCTR  techniques  based 
on  the  two  fundamental  classes  of  HRR  imaging  and  Signal  Modulation  analysis, 
developed  to  perform  target  recognition  and  identification  using  radar  systems.  A 
comparison  of  the  NCTR  techniques  and  their  relative  advantages  and  disadvantages  in 
target  recognition  performance  is  presented. 
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Chapter  III  examines  the  radar  principles  and  scattering  mechanisms  of 
electromagnetic  waves  when  a  radar  signal  intersects  any  target  in  its  propagation  path. 
A  mathematical  model  to  represent  direct  specular  and  indirect  multiple  scattering 
interactions  of  the  backscattered  radar  signal  is  presented. 

Chapter  IV  introduces  the  definitions  and  mathematical  derivations  of  the  bi¬ 
spectrum,  bi-time  and  bi-range  of  the  radar  backscattered  signal.  The  concept  of  using  the 
bi-spectral  signatures,  also  known  as  the  characteristic  bi-range  profiles  of  backscattered 
radar  signal  for  target  recognition  and  identification,  is  presented  and  studied  through  an 
example  of  a  simple  3 -points  scatterer  target. 

Chapter  V  further  investigates  the  potential  of  using  bi-spectra  technique  for  radar 
target  recognition  by  conducting  experimentation  based  on  the  Air  Force  Research 
Laboratory,  “Backhoe  Data  Dome.”  The  results  and  findings  for  three  different  test  cases 
are  presented  and  analyzed. 

Chapter  VI  discusses  the  conclusions  made  in  this  thesis  and  recommends 
objectives  for  follow-on  research  work. 
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II.  NCTR  CLASSES  AND  TECHNIQUES 


A.  FUNDAMENTAL  RADAR  TARGET  RECOGNITION  CLASSES 

NCTR  is  generally  categorized  into  two  fundamental  classes:  High  Range 
Resolution  (HRR)  imaging  and  Signal  Modulation  analysis  to  perform  radar  target 
recognition  and  identification  (Tait,  2005;  Van  Der  Heigen,  1998).  The  HRR  imaging 
class  of  NCTR  exploits  the  radar  signatures  of  a  target’s  physical  shape,  while  the  Signal 
Modulation  analysis  class  of  NCTR  exploits  the  radar  signatures  of  the  target  or  the 
dynamics  of  its  moving-parts.  A  discussion  of  the  concepts,  relative  merits  and 
limitations  of  the  current  conventional  NCTR  techniques  based  on  these  two  fundamental 
classes  is  presented  in  the  following  sections. 

B.  HRR  IMAGING 

The  basis  of  the  HRR  imaging  class  of  NCTR  is  established  on  the  ability  for  high 
radial  and  cross  range  resolutions  radar  to  resolve  the  individual  scattering  centers  of  a 
target  and  perform  radar  imaging  for  target  recognition.  NCTR  techniques  from  this 
family  include  one-dimensional  HRR  profile  or,  simply,  radial  range  profile  and  two- 
dimensional  Inverse  Synthetic  Aperture  Radar  (2D-ISAR)  imaging  (Skolnik,  2001;  Tait, 
2005;  Van  Der  Heigen,  1998). 

1.  Radial  Range  Profile 

The  radial  range  profiles  of  a  target  of  interest  essentially  represent  the  one¬ 
dimensional  images  of  the  target  as  viewed  from  different  aspect  angles.  The  dominant 
individual  scatterering  centers  or  subcomponents,  i.e.  the  parts  on  the  target  that  present 
large  Radar  Cross-Sections  (RCS)  and  subsequently  return  a  strong  radar  reflection,  are 
projected  onto  the  line-of-sight  between  the  radar  system  and  the  target. 

An  example  of  a  one-dimensional  image  that  can  be  can  be  created  by  the  HRR 
radar  system  is  illustrated  in  Figure  2.  Figure  2  shows  the  radial  range  profile  of  a  target 
generated  by  the  computer  prediction  code  using  the  Computer  Aided  Design  (CAD) 
model  of  an  aircraft  viewed  from  the  left  hand  side.  Impulse  responses,  created  using  the 
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Fourier  transform  of  the  backscattered  radar  signal  from  the  aircraft  dominant  scatterers 
(dots),  are  then  projected  onto  the  line-of-sight  in  order  to  generate  a  radial  range  profile. 


Figure  2  Radial  Range  Profile  of  a  CAD  modeled  aircraft  (Defence  Research 

and  Development  Canada). 


The  radial  range  profile  contains  information  on  the  geometry  and  heuristic 
features  of  the  target  and  can  be  used  to  infer  the  type  characteristic  of  the  target.  This 
information,  together  with  collected  reference  databases  of  measured  and  /  or  stimulated 
(i.e.  generated  using  RCS  prediction  tools  and  CAD  models)  target  signatures  of  different 
target  types  (e.g.  aircraft,  ship  and  vehicle)  and  orientations,  can  then  be  processed  by 
target  recognition  and  identification  algorithms  to  classify  the  unknown  target  into 
broader  categories,  such  as  small  jet,  large  jumbo  jet,  missile  and  tank,  etc.  down  to  the 
specific  platform  type,  such  as  F-16  fighter  jet,  Boeing-747  passenger  plane,  Kh-31  Anti- 
Radiation  Missile  and  M-l  Main  Battle  Tank,  etc. 

The  merits  and  limitations  of  the  radial  range  profile  NCTR  technique  according 
to  Van  der  Heigen  (1998)are  summarized  in  Table  1 : 


Advantages 

Disadvantages 

a.  Applicable  to  all  target  aspect  angles. 

a.  Radial  range  profile  is  highly  dependent  on 

the  target  aspect  angle.  Hence  large  data  sets  of 

target  signatures  covering  all  expected  target  types 

and  orientations  are  required. 

b.  Requires  relatively  short  time  on  target. 

Table  1  Summary  of  Radial  Range  Profile’s  Advantages  and  Disadvantages 
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2.  2D-ISAR  Imaging 

Unlike  the  radial  range  profile  which  only  uses  the  line-of-sight  or  down  range 
information  of  the  backscattered  radar  signal  for  HRR  imaging,  2D-ISAR  imaging  uses 
information  in  both  the  down  range  and  cross-range  directions  to  form  two-dimensional 
images  of  the  target  of  interests.  The  cross-range  resolution  in  2D-ISAR  can  be  thought 
of  in  terms  of  the  high  resolution  in  the  Doppler-frequency  domain.  Since  each  of  the 
major  individual  scatterer  of  a  moving  target  can  have  a  different  relative  velocity  with 
respect  to  the  radar  system,  resolution  in  Doppler-frequency,  especially  for  the  rotation 
components  of  the  target’s  motion,  will  make  it  possible  to  separate  these  scatterers  in  the 
cross-range  direction.  Figure  3  shows  an  example  of  an  ISAR  image.  The  target  is  a  B727 
jetliner  with  orientation  as  in  inset  (Borden,  1998). 


cross-rEmgcq  meters 

Figure  3  ISAR  Image  of  a  B727  Jetliner  (Borden,  1998). 
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The  merits  and  limitations  of  the  2D-ISAR  imaging  NCTR  technique  according  to 
Van  Der  Heigen  (1998)  are  summarized  in  Table  2: 


Advantages 

Disadvantages 

a.  Better  fidelity  in  mapping  the  geometry  of 

the  target  due  to  the  additional  resolution  in  the 

cross-range  dimension. 

a.  Requires  rotational  components  of  target’s 

motion  perpendicular  to  the  line-of-sight.  Hence 

target  with  weak  angular  rotation  (i.e.  pitch,  roll  and 

yaw  motions)  will  have  poor  cross-range  resolution. 

b.  More  photographic-like  2D  image  which  is 

better  suited  for  human  interpretation. 

b.  Requires  complex  motion  compensation 

and  auto-focusing  algorithm  of  2D-ISAR  images 

which  is  computationally  heavy. 

c.  Limited  knowledge  on  the  precise  project 

plane  for  the  rotation  vectors,  making  recognition 

complicated. 

Table  2  Summary  of  2D-ISAR  Imaging’s  Advantages  and  Disadvantages 


C.  SIGNAL  MODULATION  ANALYSIS 

The  Signal  Modulation  analysis  class  of  NCTR  is  based  on  the  characterizations 
of  the  reflected  radar  modulation  signals  from  the  rotating  parts  of  the  engine,  such  as  the 
compressor  blades  in  a  jet  engine,  propellers  on  a  propeller-aircraft  and  rotors  on  a 
helicopter,  in  the  Doppler  domain.  The  NCTR  techniques  from  this  family  are 
prominently  used  to  perform  target  recognition  and  the  identification  of  aircrafts  and 
include  the  following:  Jet  Engine  Modulation  (JEM),  Propeller  ROtor  Modulation 
(PROM)  and  HElicopter  Rotor  Modulation  (HERM)  (Skolnik,  2001;  Tait,  2005;  Van  Der 
Heigen,  1998). 


1.  JEM,  PROM  and  HERM 

The  compressor  blades  in  a  jet  engine,  aircraft  propellers  and  rotating  blades  and 
hub  of  a  helicopter,  each  produces  a  unique  and  distinctive  modulation  of  the  radar  target 
signatures.  These  modulations  in  the  characteristic  spectrums  of  JEM,  PROM  and 
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HERM  can  be  analysed,  and  attributes  pertaining  to  the  information  on  engines  / 
propellers  /  rotors  types,  quantities,  number  of  blades,  ranges  of  rotation  frequencies,  etc, 
can  be  extracted.  Examples  of  the  JEM  spectrum  of  a  twin-engine  aircraft,  PROM 
spectrum  of  a  DC-7  four-engine  aircraft  and  HERM  spectrum  of  a  hovering  helicopter 
are  shown  in  Figures  4,  5  and  6  respectively. 


Figure  4  JEM  Spectrum  of  a  Twin-Engine  Aircraft  (Van  Der  Heigen,  1998) 


Doppler  frequency  (Hz) 

Figure  5  PROM  Spectrum  of  a  DC-7  four-engine  aircraft  (Skolnik,  2001) 
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Figure  6  HERM  of  a  Hovering  Helicopter  (Bullard  &  Dowdy,  1991) 


Since  each  aircraft  has  a  unique  engine,  propeller  or  rotor,  these  characteristic 
parameters  can  be  matched  against  a  reference  database  or  table  that  contains 
specifications  for  a  variety  of  target  engines  /  propellers  /  rotor  types  for  engines 
associated  target  recognition  and  identification. 

The  merits  and  limitations  of  JEM,  PROM  and  HERM  analysis  NCTR  techniques 
according  to  Van  Der  Heigen  (1998)  are  summarized  in  Table  3: 


Advantages 

Disadvantages 

a.  Small  and  simple  target  database  since  it  is 

comprised  of  engine  specification  rather  than  radar 

target  signature  or  image. 

a.  Dependent  on  the  target  aspect  angle. 

Useful  spectrums  mostly  in  the  head-on  or  tail-on 

orientation. 

b.  Requires  relatively  short  time  on  target. 

b.  Requires  strong  signal-to-noise  ratio  for 

effective  classification.  Hence  techniques  limited  to 

relatively  short  ranges 

c.  Fast  classification  since  almost  all  aircrafts 

have  a  unique  engine  type. 

Table  3  Summary  of  JEM,  PROM  and  HERM  Analysis’s  Advantages  and  Disadvantages 
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III.  RADAR  PRINCIPLES  AND  TARGET  SCATTERING 


A.  OVERVIEW  OF  RADAR  CONCEPT 

RAdio  Detection  And  Ranging  (RADAR)  is  a  method  of  using  electromagnetic 
waves  to  remotely  sense  the  presence  and  extract  information  pertaining  to  the  position, 
velocity  and  identifying  characteristics  of  targets  of  interest.  An  overview  of  the  basic 
radar  operating  principle  and  electromagnetic  wave  propagation  theory  is  presented  in  the 
following  sections. 

1.  Principle  of  Radar  Operation 

The  basic  radar  block  diagram  and  principle  of  radar  operation  is  illustrated  in 


transmitter duplexer antenna  •-  electromagnetic  wave-^aim 

display  receiver^- duplexer  —  antenna  —  echo  signal  —  aim 

Figure  7  Basic  Radar  Block  Diagram  and  Operating  Principle  (Wolff,  1997) 

Using  a  transmitter,  the  radar  system  first  generates  an  electromagnetic  signal 
(such  as  a  short  sinusoidal  pulse)  that  is  radiated  into  free  space  through  a  duplexer  and 
an  antenna.  The  duplexer  switches  the  antenna  of  monostatic  single-antenna  systems 
between  the  transmitter  and  receiver,  allowing  the  antenna  to  be  shared  for  both 
transmitting  and  receiving  functions,  while  the  antenna  is  used  to  concentrate  and  direct 
the  electromagnetic  energy  towards  the  target. 
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As  the  electromagnetic  wave  propagates  in  free  space,  a  fraction  of  the 
transmitted  electromagnetic  energy  is  intercepted  by  the  target  and  scattered  in  all 
directions.  The  reflected  energy  of  the  echo  signal  directed  back  towards  the  radar 
system,  which  is  also  known  as  the  radar  backscattered  signal,  is  then  collected  by  the 
same  antenna  and  redirected  to  the  receiver  through  the  duplexer.  Thereafter,  the 
receiver  performs  the  signal  processing  of  the  received  echo  signal  in  order  to  detect  the 
presence  of  the  target  and  extract  information  pertaining  to  its  position,  velocity  and 
identifying  characteristics.  The  process  is  completed  with  the  presentation  of  this 
information  on  an  operator  display. 

2.  Theory  of  Electromagnetic  Wave  Propagation 

The  fundamental  laws  of  electricity  and  magnetism,  also  known  as  Maxwell’s 
equations,  predicted  the  phenomenon  of  electromagnetic  waves  in  the  form  of  sinusoidal 
electric  E-field  and  magnetic  H-field,  which  are  oscillating  at  the  wave’s  frequency  and 
mutually  orthogonal  to  the  direction  of  propagation,  as  depicted  in  Figure  8. 

Direction  of 


Figure  8  Electromagnetic  Wave  With  Mutually  Orthogonal  E-Field,  H-Field 

and  Direction  of  Propagation  Components 


16 


By  convention,  the  electric  E-field  is  usually  adopted  in  the  representation  of  the 
electromagnetic  wave.  A  derivation  of  the  electric  E-field  wave  equation  from 
Maxwell’s  equations  is  presented  according  to  Sadiku’s  Elements  of  Electromagnetics 
(2001).  Beginning  from  Maxwell’s  equations  for  a  vacuum: 


V-E  =  0 

dt 

V  H  =  0 

d 

VxH  =  £,,  —  E 

“at 


where 

V  is  the  vector  differential  operator, 

jUo  =  4n  x  10'7  H/m  is  the  penneability  of  the  vacuum  and 
12 

So  =  8.854  x  10'  ‘  F/m  is  the  permittivity  of  the  vacuum. 


(3.1) 

(3.2) 

(3.3) 

(3.4) 


Taking  the  curl  of  equation  (3.2): 


V  x(Vxf)  =  V : 


-«4" 


(3.5) 


Evaluating  the  left  hand  side  of  equation  (3.5)  with  the  vector  identity  of 
Vx(Vxd)  =  V(V-d)- V2A  and  simplifying  the  expression  using  equation  (3.1), 


Vx(Vx£)  =  V(V-E)-V2E  =  -V2E 


(3.6) 


Evaluating  the  right  hand  side  of  equation  (3.5)  and  simplifying  the  expression  using 
equation  (3.4), 


Vx 


-«4 H 


d 


d2 


=  -Mo  ^(V  x  H)  =  -MAo 


(3.7) 
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Equating  equations  (3.6)  and  (3.7),  results  in  a  differential  equation  for  the  electric  E- 
field  wave  equation, 

2  52 

V2E-<u0£0  —  E  =  Q  (3.8) 

dt2 

where 

fioSo  =  Me"  and 

c  =  2.998  x  108  m/sec  is  the  speed  of  light. 

3.  Radar  Signal 

By  solving  the  electric  E-field  wave  equation  in  equation  (3.8),  the  radar  signal 
for  an  electromagnetic  wave  propagating  from  a  radar  system  can  be  represented  as  a 
general  fonn  of  a  complex  sinusoid  electric  E-field  varying  in  both  time  and  space  and 
expressed  as 

E(r,t)  =  a(oj)[cos(ojI  -k-r)  —  i  sin(cot  -  k  ■  r)]  (3.9) 

where 

A  is  the  amplitude  of  the  E-field  signal, 

co  =27rfis  the  radial  frequency, 

k  =  —  is  the  wave  number, 

A 

A  is  the  wavelength, 

r  is  the  range  from  the  radar  system  and 
t  is  the  time  duration  elapsed. 

Elsing  Euler’s  fonnula,  e~lx  =  cosx-i  sinx ,  the  radar  signal  at  a  given  position  and  time, 
E(r,  t )  can  be  rewritten  in  the  exponential  fonn  as 

E(r,t)=A(co)e^cot~k'r^  (3.10) 
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To  compute  the  echo  signal  collected  at  the  radar  receiver,  the  range  from  the  radar 
system  is  set  to  zero,  r  =  0,  hence  the  received  radar  signal  in  the  time  domain  is 

E(t)  =  A(co)e~i^  (3.11) 

While  the  corresponding  received  radar  signal  in  the  frequency  domain  is 

H(<d)=  A(co)e~i^  (3.12) 


B.  TARGET  SCATTERING  MECHANISMS 

The  mechanisms  by  which  target  scattering  occurs,  are  important  considerations 
for  determining  the  physical  size,  shape,  orientation  and  material  composition  of  the 
target,  which  then  account  for  the  portion  of  electromagnetic  energy  that  is  reflected  back 
in  the  direction  of  the  radar  system  and  consequently  the  probability  of  target  detection, 
recognition  and  identification. 

As  discussed  in  Tait  (2005),  the  main  interactions  of  the  electromagnetic  wave 
with  the  target  are:  reflection;  absorption;  refraction;  diffraction  and  creeping  surface 
waves.  Figure  9  illustrates  the  interaction  effects  of  reflection,  absorption  and  refraction. 
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At  the  surface  boundary  of  the  target,  the  incident  electromagnetic  energy  is 
reflected,  refracted  and  absorbed  by  the  interactions  of  the  radar  wave  with  the  target.  It 
can  be  observed  that  the  fraction  of  incident  electromagnetic  energy  reflected  back 
towards  the  radar  system  is  dependent  on  the  portion  of  energy  that  is  absorbed  and 
refracted  into  the  material,  as  well  as  the  angle  of  incidence  to  the  material  surface. 

Another  fonn  of  scattering  mechanism,  through  the  interaction  of  the  radar  wave 
with  the  target,  is  diffraction.  This  knife-edge  effect  is  explained  by  Huygens'  Principle, 
which  states  that  a  well-defined  obstruction  to  an  electromagnetic  wave  acts  as  a 
secondary  source  and  creates  a  new  wavefront.  This  new  wavefront  then  propagates  into 
the  geometric  optics  shadow  area  of  the  obstacle,  enabling  the  radar  wave  to  bend  around 
the  edges  of  the  target,  propagate  into  the  target’s  structure  through  openings  and  detect 
the  interior  parts  of  the  target  that  would  be  considered  as  non  line-of-sight  as  shown  in 
Figure  10. 


Figure  10  Huygens’  Principle  of  Diffraction  (Potter) 
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Finally  the  creeping  surface  wave  interaction  effect  allows  for  the  propagation  of 
radar  electromagnetic  wave  along  the  surface  of  a  target  using  surface  wave  to  access  the 
geometric  optics’  blind  region. 

C.  REPRESENTATION  OF  RADAR  BACKSCATTERED  SIGNAL 

Targets  are  generally  composed  of  many  individual  discrete  specular  scatterers, 
each  with  different  scattering  properties  which  correspond  to  a  physical  feature  of  the 
target,  for  example  the  nose,  engines,  main  and  tail  wings  of  an  aircraft. 

Each  scatterer,  when  intercepted  by  a  radar  electromagnetic  wave,  will  reflect  an 
echo  signal  characterized  by  an  amplitude  and  a  phase  back  towards  the  direction  of  the 
radar  system,  which  will  then  sum  up  to  a  resultant  radar  backscattered  signal  at  the  radar 
receiver.  In  addition,  there  will  be  contributions  from  the  interactions  between  the 
multiple  individual  scatterers  that  will  further  complicate  the  target  characteristics. 

A  simple  mathematical  model  for  the  direct  specular  and  indirect  multiple 
scattering  interactions  of  the  backscattered  radar  signal  is  discussed  and  presented  in  the 
following  sections  as  referenced  from  Borden  (1998)  and  Jouny  (1994). 

1.  Mathematical  Model  for  Direct  Specular  Scattering 

The  radar  backscattered  signal  is  assumed  to  be  composed  of  two  major 
components.  The  first  component  represents  the  specular  direct  scattering  response  from 
the  individual  scatterers,  distributed  over  the  target.  The  propagation  path  of  the  signal 
can  be  represented  in  Figure  1 1  as: 


tp/2  one-way  propagation  time 


pth  scatter  point 


tp/2  one-way  propagation  time 


Radar 


Round-trip  propagation  time  =  tp/2  +  tp/2  =  tp 


Figure  1 1  Propagation  Path  of  Direct  Specular  Scattering  Response  from 

Individual  Scatterers 
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Using  equation  3.12,  this  component  of  the  backscatter  signal  can  be  expressed 
mathematically  as  a  summation  of  the  reflected  echo  signals  from  all  of  the  individual 
discrete  specular  scatterers: 


where 


hM=  ±Ap(ay‘^°*') 

P= 1 


(3.13) 


m  is  the  number  of  scattering  centers, 

Ap  is  the  amplitude  of  the  scattered  signal  from  the  /;th  scattering  centre, 

y/p  is  the  phase  shift  due  to  pi\\  scatterer  and  statistically  characterized  by  the 
uniform  probability  density  distribution  over  the  interval  (0,  2ti), 

co  is  the  radial  frequency  (=2/z/)  and 

tp  is  the  round-trip  propagation  time  between  the  radar  and  the  /till  scattering 
centre. 


2.  Mathematical  Model  for  Indirect  Multiple  Scatter  Interactions 

The  second  major  component  of  the  radar  backscattered  signal  represents  the 
effect  of  interaction  between  multiple  scattering  points.  The  1st  order  interaction  between 
any  two  scattering  points  is  represented  in  Figure  12  as: 


& 

Radar 


tp/2  one-way  propagation  time 


tpq  interaction  time 


pth  scatter  point 


tq/2  one-way  propagation  time 


qth  scatter  point 


Round-trip  propagation  time  =  tp/2+(tpq)+tq/2=  (tp+tq)/2  +  tpq 

Figure  12  Propagation  Path  of  Indirect  Scattering  Response  from  Multiple 

Scatterers’  Interactions 
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Using  equation  3.12,  this  second  component  can  also  possibly  be  modeled 
mathematically  as  a  summation  of  the  reflected  echo  signals  from  all  of  the  multiple 
scatterers’  interactions: 

m  m  ~l 

YJApq((°)e 

p= 1  q=\q*p 

+  ...higher  -  order  interaction 

where 

Apq  is  the  amount  of  coupling  between  the  p\h  and  <yth  scattering  points, 

y/q  is  the  phase  (uniformly  distributed  over  2 71)  of  the  of  c/th  scatter, 

tq  is  the  round-trip  propagation  time  between  the  radar  and  the  c/th  scattering 
centers  and 

tPq  is  the  interaction  time  between  the  /;th  and  c/th  scattering  centers. 

Since  the  interaction  time,  tpq  between  the  /cth  and  c/th  scattering  centers  is 
generally  unknown  for  complex  targets,  it  is  assumed  to  be  equivalent  to  (tp+tq)/2  for  the 
present  (although  this  can  be  easily  generalized).  The  radar  backscattered  signal  due  to 
multiple  interactions  is  then  scattered  at  a  round-trip  propagation  time  of  (tp+tq)/2+(tp+tq)/2 
=  tp+tq  and  the  second  major  component  of  this  signal  can  be  rewritten  as: 


v+v+w 


t 


L+t 


mm  .r  /  \\ 

P= 1  q=\,q*p 

+ . . . .  higher  -  order  interaction 


(3.15) 


3.  Mathematical  Model  for  Combined  Direct  Specular  and  Indirect 
Multiple  Scattering  Interactions 

Collectively,  the  coherent  radar  backscattered  signal  of  the  target,  H(a)  at  radial 
frequency,  (o=27jf,  is  a  combination  of  both  the  direct  specular  and  multiple  scatter 
interaction  terms  of  Hj(a>)  and  H2(<x>). 

Ignoring  the  effect  of  higher-order  interactions  and  assuming  that  the  signal 

amplitudes  are  frequency-independent,  the  radar  backscattered  signal,  H(a)  is  given  as: 
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(3.16) 


H(a>) 


=  Hl(a>)  +  H2(a>) 


1 4 

P= 1 


-i^+a*,) 


+ 


z  z^ 

p=lq=l,q±p 


-r’[v/+v/+°j(tl+tq)] 


pq 


However  in  target  range  profiling,  the  point  scatterer  data  are  often  observed  in 
terms  of  relative  range  instead  of  the  round-trip  propagation  times.  As  such,  the  radar 
backscattered  signal  can  be  rewritten  using  the  radar  range  relationship  of  t=2r/c  as: 


Htf) 


m  - 

\A7tR—+yf\  m 

m  —i 

^^{RP+Rq)—+¥P+Yq 

Zv 

l  c  J  +  £ 

lLi^pqe 

c 

p= i 

p=\q- 

=l,q*p 

(3.17) 


where 


Rp  is  the  range  between  the  radar  and  the  pth  scatter  point, 

Rq  is  the  range  between  the  radar  and  the  <yth  scatter  point, 

/  is  the  frequency  and 

c  is  the  speed  of  light. 


D.  EXAMPLE:  RADAR  BACKSCATTERED  SIGNAL  FROM  A  SIMPLE  3- 
POINT  SCATTERERS  TARGET 

To  illustrate  mathematical  modeling  of  the  backscattered  radar  signal,  consider  a 
3-point  scatterers  hypothetical  target  with  scatter  points  at  ranges  of  Ri,  Rj  and  Rj  from 
the  radar  as  shown  in  Figure  13. 


Range,  R3 


Figure  13  Direct  and  Indirect  Propagation  Path  of  a  3-Point  Scatterers  Target 


24 


Using  equation  3.17,  the  radar  backscatter  signal  for  a  given  frequency, /is  given 


as: 


H(f)  =  2U 


3  -/  4kR „—+>/„  I  3  3 


X  I  -V 


4n{Rp+Rq)UVp+Vq 


P= 1 


/e 


/ 


+412e 


+  42e 

f 

-/I  4;r(.fi1  +,ff2 )— +^]  +i//2 


p= i  q=i,q*p 
-if  4nR2— +y/2 


+ 


+zf, 


r 

-il  4a-(S[+i?2)— +V/i+!<,2 


+4/ 


f 

— /I  - 


+  ^3^ 


^23^ 


^32  ^ 


f 

-il  47iRi  — +(/3 
c 

f 

4  rr(Rl  +R2  )—+y/ \  +i//2 

f 

4x(R2+R}  )l-+y2+<//2 
j^43-(i?2+i?3  )^-+(/2+!^3 


When  R[  =  0  m,  7/  =  1.5  m  and  //  =  3.6  m,  the  corresponding  radar  backscatter 
signal  can  be  rewritten  as: 


H(f)  =  C,  +  C2e 


i\  4^(1. 5)— 

—i 

4^(3.6)C 

-(4^5.1)^ 

L 

+  C-,6 

c  _ 

+  CAe 

where 


Ci,  C2,  C3  and  C4  are  constants. 


E.  EXAMPLE:  RADIAL  RANGE  PROFILE  OF  A  SIMPLE  3-POINT 

SCATTERERS  TARGET 

To  better  understand  and  visualize  the  contributions  of  individual  scatterers  and 
their  multiple  interactions  in  the  radar  backscattered  signal,  the  radial  range  profile  of  the 
target  is  examined.  The  radial  range  profile  of  any  target  can  be  generated  by  evaluating 
the  amplitudes  and  phases  of  the  radar  backscattered  signals  over  a  frequency  band  of 
evenly  spaced  frequency  intervals  and  performing  Fourier  Transform  to  convert  the 
signals  into  the  time  domain  (Jouny,  1994).  The  radial  range  profile  is  a  plot  of  the 
magnitude  of  the  impulse  response  versus  the  projected  down  range  along  the  line-of- 
sight  between  the  radar  system  and  the  target,  where  scatterings  from  the  individual 
specular  scatterers  and  their  multiple  interactions  will  appear  as  peak  responses. 
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The  radial  range  profile  of  the  simple  3-point  scatterers  target,  computed  using 
MATLAB’s  Fast  Fourier  Transfonn  of  the  associated  radar  backscattered  signals  over  a 
frequency  band  from  5.0  GHz  to  6.0  GHz,  is  illustrated  in  Figure  14.  The  related 
MATLAB  program,  “backscatter  signal.m”  and  source  code  used  to  generate  this  plot 
can  be  found  in  Appendix  A. 


Figure  14  Radial  Range  Profile  of  3-Point  Scatterers  Target 


From  Figure  14,  it  can  be  observed  that  the  first  three  peak  impulse  responses  are 
consistent  with  the  relative  positions  and  geometry  of  the  3-points  scatterers  at  ranges  Ri 
=  0  m,  R2  =  1.5  m  and  R3  =  3.6  m,  while  the  fourth  peak  impulse  response  is  a  “ghost” 
artifact  that  does  not  map  to  any  scatterers  and  represents  the  collective  summation  of  the 
multiple  scatters  interaction  terms  between  the  2nd  and  3ld  scattering  points  at  a  range 
equivalent  to  R2  +  fG  =  1.5  m  +  3.6  m  =  5.1  m.  (Borden,  1998;  Jouny,  1994). 
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IV.  BI-SPECTRAL  NCTR  TECHNIQUE 


A.  OVERVIEW  OF  BI-SPECTAL  METHOD 

While  the  radial  range  profile  NCTR  technique  proves  to  be  a  viable  means  of 
performing  target  recognition,  its  application  is  predominantly  based  on  the  correct 
mapping  and  correlation  of  the  peak  impulse  responses  from  the  direct  specular  scatterers 
to  the  actual  physical  layout  of  the  target  structure  and  its  characteristic  features. 

However,  as  discussed  and  observed  in  the  last  section  of  Chapter  III,  some  of 
these  peak  impulse  responses  in  the  radial  range  profile  do  not  correspond  to  any  target 
structures  and  represent  the  multiple  interactions  between  major  individual  scatterers 
instead.  Thus,  a  bi-spectral  method,  to  differentiate  the  direct  specular  scattering 
responses  from  the  multiple  scatterers  interaction  responses,  is  proposed  in  Jouny  (1994), 
Botha  &  Spoelstra  (1992),  Garber,  Jouny  &  Moses  (1995)  and  Jounty  &  Walton  (1990) 
as  another  HRR  imaging  class  of  NCTR  technique  to  augment  the  radial  range  profile  in 
performing  radar  target  recognition  and  identification. 

1.  Definition  of  Bi-Spectrum 

In  Higher-Orders  Statistic  (HOS)  signal  processing,  the  bi-spectrum  of  a 
continuous-time  signal,  h(t)  is  defined  as  the  Fourier  Transform  of  the  third-order 
cumulant-generating  function  and  is  given  in  Jouny  (1994),  Botha  &  Spoelstra  (1992), 
Garber,  Jouny  &  Moses  (1995)  and  Jounty  &  Walton  (1990): 

GO  OO 

B(col,co2)  =  j  J  C3  {t\,t2  ^dt\dt2  C4-1) 

—  00  —  00 

where 

C3  represents  the  triple  correlation  of  the  signal  and  is  expressed  as: 

C2>(tl,t2)  =  J \h*(t)h(t  +  tl)h(t  +  t2)dt  =  e\i*  (t)h(t  +  tx)h(t  +  t2)\  (4-2) 

—  00 

E{  }  denoted  the  expectation  and 

*  superscript  asterisk  denoted  a  complex  conjugate. 
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2.  Definition  of  Bi-Time 

However,  in  radar  signal  processing,  the  radar  backscattered  signals  are  often 
processed  and  analyzed  as  discrete  signals  in  the  frequency-domain  (i.e.  expressed  as  a 
function  of  frequency).  The  third-order  cumulant  or  triple  correlation  of  the  frequency- 
domain  signal,  H(a>)  is  then  defined  by  Jouny  (1994),  Botha  &  Spoelstra  (1992),  Garber, 
Jouny  &  Moses  (1995)  and  Jounty  &  Walton  (1990)  as: 

00 

C3  {(O\,C0  2  )  =  |  H*  (g>)h(cO  +  (0\  )H(C0  +  0l>2  )dco 

—  00  (4-3) 

+  C0\)H(C0  +  fi>2)} 

The  equivalent  bi-spectrum,  expressed  as  a  function  of  time  after  applying  the 
Fourier  Transform,  is  known  as  the  bi-time.  The  bi-time  of  a  discrete  time-domain  bi- 
spectral  signature  is  then  expressed  as: 

00  00  ,  . 

X  XC3(^,fi^>_,W‘+<W  (4.4) 

co{  —  00  co2  —  00 

and  can  be  represented  in  terms  of  impulse  response  as: 

B(t],t2)=  )h{ti )h* (t]  +t2)\  (4.5) 

where 

h( t)  is  the  impulse  response  as  a  function  of  time 


3.  Definition  of  Bi-Range 

The  bi-time  can  be  further  expressed  as  a  function  of  range  using  the  radar  range 
and  round-trip  propagation  time  relation  of: 


r  = 


ct 

2" 


t 


(4.6) 


where 

r  is  the  range  from  the  radar  to  the  target, 

t  is  the  round-trip  propagation  time  needed  for  the  radar  signal  to  propagate 

to  the  target  and  back,  and 
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c 


is  the  speed  of  light. 


Substituting  t=2r/c  and  /  =o.)/2n  into  the  expression  for  bi-time,  derives  the 
corresponding  bi-range  as  (Jouny  (1994),  Botha  &  Spoelstra  (1992),  Garber,  Jouny  & 
Moses  (1995)  and  Jounty  &  Walton  (1990): 


00  00 


B{n  ,n)  =  X  Xc3(/b/2> 

f\  =  ~ 00  f  2  =  ~°° 


f  2r.  )  1 

(  2  r2  V 

2< 

—  +2#J 

l  c  )  1 

V  c  )_ 

00  00 


=  S 

/l=-CO/2=-CO 


C/ih+M) 

C 


(4.7) 


and  bi-range  can  then  be  represented  in  terms  of  impulse  response  as: 

£(/1  ,g  )=  e \ll(r\  )^(r2 K (n  +G)}  (4-8) 

where 


/z(r)  is  the  impulse  response  as  a  function  of  range. 

From  equation  4.7,  it  can  be  deduced  that  if  the  range  profde  for  a  given  radar 
target  contains  three  responses,  h(r\),  h{ri)  and  h(r 3),  two  of  which,  h(r\)  and  h(ri)  are 
due  to  direct  specular  scattering,  and  one  is  due  to  an  interaction  between  the  direct 
scatterers,  h{ri)  ,  the  bi-range  will  only  detect  this  tenn  if  the  interaction  tenn  appears  at  a 
distance  from  the  radar  equal  to  the  sum  of  the  distances  of  the  direct  terms,  i.e.  at  rj,  =  r\ 
+  ro  according  to  Jouny  (1994). 


B.  EXAMPLE:  BI-RANGE  OF  A  SIMPLE  3-POINT  SCATTERERS  TARGET 

Using  the  same  simple  3-point  scatterers  target  example  as  described  in  Chapter 
III,  the  bi-range  of  the  target  can  be  computed  mathematically  using  equation  4.7  as: 

B(n  ,r2)  =  E  )h(r2  )h  *  (r{  +  r2  )} 

44  r  /  \  [  vi  (4.9) 

=  X  X  ^P^q^pq  —  Bp  5  r2  —  Bq  j"1"  l  —  Bq  ’  r2  —  Bp  ). 

p=\q=\ 
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Expanding  the  tenns  for  the  above  expression: 


B(r\ > r2 )  -  AAA  1  [^(rl  -  A  , r2  ~  A  )  +  S(rx  -Rl,r2-Rl )] 

+  A\A2  A\2  \s(rx  ~R\,r2  -R2)  +  S{r\  -  R2 ,  r2  -  Rx )] 

+  A1A3A13  \s{r\  -Ri,r2-R3)+  S(rx  -R3,r2-Rx )] 

+  A\  A4  A\ 4 \8{rx  -R\,r2  - R4 )  +  8(r\  - R4 , r2  -  R\ )] 

+  A2  A\ A2X \S(rx  -R2,r2  -R\)  +  S{j\  -  Rx , r2  -  R2 )] 

+  A2  A2  A22  [c)(n  - R2,r2  - R2)+  S{i\  -  R2,r2  -  R2 )] 

+  A2A3A23[S(ri  -  R2 , r2  -R3)  +  S(r\  -R3,r2  -  R2)\ 

+  A2A4A24[A(ri  -R2,r2  -  R4)  +  d(r\  -R4,r2  -  R2)\ 

+  A3AlA3l[d'(rl  -  R3 ,  r2  -  R\  )  +  S{r]  -R\,r2  -  R3 )] 

+  A3A2A32[S(rl  -  R3,r2  -R2)+8{rx  -R2,r2  -^3)] 

+  AAA3\d(rl  ~R3’r2  -  R3  )  +  Ar\  -  A  > r2  -^3)] 

+  A3 A4 A34  [S{r\  -  R3,r2  - R4)+  S(rx  -  R4 , r2  -  R3 )] 

+  A4 A\ A4 1  \s{rx  -R4,r2  ~R\)  +  S{i\  -  Rx , r2  -  R4 )] 

+  A4  A2  A42  \S{rx  -  R4,r2  -  R2)+  S(rx  -  R2 , r2  -  R4 )] 

+  A4A3A43\s(r\  -  R4,r2  -R3)+S(rx  ~R3,r2  -^4)] 

+  A4 A4 A44 [s(rx  - R4,r2  - R4)  +  S{rx  - R4,r2  - R4 )] 

Substituting  Ri  =  0  m,  R2  =1.5  m,  R<  =  3.6  m  and  R4  =5.1  m,  the  bi-range  can 
be  rewritten  as: 

S(ri.r2)  =  C1^'-2)]+C2,[<S(r1  -1.5, r2  -1.5)] 

+  C3  ’  [d(rx  -  3 .6,  r2  -  3 .6)]  +  C4  ’  [S(rx  -  5. 1,  r2  -  5.  l)] 

+  C5  ’  [4n ,r2 -\. 5)  +  d{rx  -1.5,  r2 )] 

+  C6  ’  [S(rx  ,r2-  3.6)  +  S{rx  -  3 .6,  r2 )] 

+  C7'[s(rl,r2-5.i)  +  S(rl-5A,r2)\ 

+  C8  ’  [S(rx  - 1 .5,  r2  -  3 .6)  +  S(r\  -  3 .6,  r2  - 1 .5)] 

+  Cq'\s{ri  -1.5 ,r2  —  5 . l)  +  -5.1,r2  -1.5)] 

+  C]0'[s(r]  -3.6 ,r2  -5.l)  +  fi(/i  -5.1, r2  -3.6)] 

where  C{,  C2\  C3\  C4\  C5\  C6\  C7\  C%\  C9’  and  CV  are  constants. 

Graphically,  the  bi-range  of  the  3 -point  scatterers  target  can  be  computed  using 

MATLAB’s  two-dimensional  (2-D)  discrete  Fourier  Transform  of  the  third-order 

cumulant  function  of  the  radar  backscatter  signal.  The  results  are  then  plotted  as  the 

three-dimensional  (3-D)  and  2-D  (contour  plot)  bi-range  profiles  of  the  3 -point  scatterers 
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target  in  Figure  15  and  16  respectively.  The  related  MATLAB  programs,  “bispectral.m” 
and  source  code  used  to  generate  these  plots  can  be  found  in  Appendix  B. 


Bi- range,  r2  (m)  Bi-range,  r,  (m) 

Figure  15  3-D  Bi-Range  Profile  of  3-Point  Scatterers  Target 


Figure  16  2-D  Bi-Range  Profile  of  3-Point  Scatterers  Target 
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From  Figures  15  and  16,  it  can  be  observed  that  the  bi-range  profiles  are 
generated  by  responses  at  the  range  pairs  (Rp  and  Rq)  for  p  and  q=  1,. ...  4  with  amplitude 
of  ApAqApq.  The  locations  of  the  responses  are  also  noted  to  be  consistent  with  the 
theoretically  computed  positions. 

C.  CHARACTERISTIC  BI-RANGE  PROFILE 

Since  a  non-zero  bi-range  response  is  an  outcome  arising  from  the  correlated 
interaction  between  the  responses  at  range  r\  and  ri,  it  is  clear  that  the  peak  responses  in 
the  bi-range  profile  of  the  radar  backscattered  signals  represent  strong  additive  multiple 
interactions  between  the  scatterers  or  subcomponents  of  the  target. 

Therefore,  the  characteristic  bi-range  profile  of  a  target,  which  highlights  the 
defining  peak  responses  (within  an  order  of  magnitude)  at  coordinates  (r\,  1-2)  can 
potentially  be  employed  to  characterize  the  unique  composition  and  layout  of  the  target 
for  recognition  and  identification.  The  characteristic  bi-range  profile,  showcasing  the 
peak  responses  of  the  3-points-scatterer  target  example  is  illustrated  in  Figure  17. 


Figure  17  2-D  Characteristic  Bi-Range  Profile  of  3-Point  Scatterers  Target 
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It  can  be  observed  from  Figure  17  that  the  peak  responses  in  the  bi-range  profile 
of  the  3-points  scatterer  target  at  (r\,  ri)  =  (1.5m,  3.6m)  corresponded  to  the  resulting 
multiple  interactions  between  the  2nd  scattering  point  at  1.5m  and  3ld  scattering  point  at 
3.6m,  as  expected. 

The  merits  and  limitations  of  bi-range  profile  NCTR  technique  according  to  Bao, 
Shi  &  Zhang  (2001)  and  Bao,  Pei  &  Xing  (2001)  are  summarized  in  Table  4: 


Advantages 

Disadvantages 

a.  Retains  both  the  amplitude  and  phase 

information  of  a  signal.  Hence  bi-spectral 

signatures  are  translation  (or  shift)  invariant  and 

robust  to  aspect  dependence. 

a.  Requires  intensive  computation  of  two- 

dimensional  high  range  resolution  data.  The 

b.  Good  clutter  and  noise  rejection 

performance  due  to  ability  to  statistically  suppress 

the  clutter  or  noise  characterized  by  zero-mean 

Gaussian  probability  density  distribution. 

number  of  bi-spectral  data  is  the  square  of  the  HRR 

profile. 

Table  4  Summary  of  Bi-range  Profi 

Ic’s  Advantages  and  Disadvantages 
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V.  BI-SPECTRAL  NCTR  EXPERIMENTATION 


A.  “BACKHOE  DATA  DOME”  EXPERIMENTATION 

To  further  examine  the  potential  of  using  the  characteristic  bi-range  profile  of  a 
complex  target  for  NCTR  and  identification,  experimentation  based  on  the  Air  Force 
Research  Laboratory  (AFRL),  “Backhoe  Data  Dome,  Version  1.0”  by  Lin  &  Naidu 
(2004)  are  performed  and  presented  in  the  following  sections. 

1.  AFRL  “Backhoe  Data  Dome  Version  1.0” 

The  “Backhoe  Data  Dome,  Version  1.0”  consists  of  simulated  wideband  (7-13 
GHz),  full  polarization,  complex  backscatter  data  from  a  backhoe  loader  generated  using 
a  Computer- Aided-Design  (CAD)  model  in  free  space  as  shown  in  Figure  18.  The 
MATLAB  program  “plotfacet.m”  available  from  the  Air  Force  Research  Laboratory 
Sensor  Data  Management  System  was  used  to  plot  all  3-D  graphics  of  the  backhoe 
loader. 


Figure  18  CAD  Model  of  Backhoe  Loader 


The  backscatter  data  are  then  synthetically  generated  over  a  full  2n  steradian 
viewing  hemisphere  with  14  samples  per  degree  in  both  azimuth  and  elevation.  There  is 
one  sample  every  11.75  MHz  over  the  6  GHz  bandwidth  and  full  polarizations;  i.e. 
horizontal  linear  polarization  on  transmit  and  receive  (HH),  vertical  linear  polarization  on 
transmit  and  receive  (VV)  and  horizontal  linear  polarization  on  transmit  and  vertical 
linear  polarization  on  receive  (HV). 
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A  representation  of  the  k-space  data  dome  is  as  shown  in  Figure  19. 
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Figure  19  Backhoe  Data  Dome  Representation  in  k-Space  (Lin  &  Naidu,  2004) 

In  this  experiment,  two  sets  of  2-D  backhoe  backscattered  data  which  contained 
VV,  HH  and  HV  polarized  k-space  data  from  a  110  degree  azimuth  cut  (from  350 
degrees  to  100  degrees  in  azimuth)  available  through  Lin  &  Naidu  (2004)  are  used,  one  at 
0  degree  elevation  and  the  other  at  30  degree  elevation. 

With  the  method  described  in  the  preceding  Chapters,  characteristic  bi-range 
profiles  analysis  on  the  2-D  k-space  data  is  perfonned  for  three  different  test  cases,  with 
the  backhoe  loader  setup  in  the  following  orientation  and  using  the  specified  polarization: 

a.  Test  Case  1:  Orientation  of  backhoe  loader  at  0  degree  azimuth  and  0 
degree  elevation  with  W,  HH  and  HV  polarizations. 

b.  Test  Case  2:  Orientation  of  backhoe  loader  at  90  degree  azimuth  and  0 
degree  elevation  with  VV,  HH  and  HV  polarizations. 

c.  Test  Case  3:  Orientation  of  backhoe  loader  at  60  degree  azimuth  and  30 
degree  elevation  with  VV,  HH  and  HV  polarizations. 
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2.  Test  Case  1:  Setup  and  Results  (0azimuth=O°,  (peievation=0o) 

In  Test  case  1,  the  backhoe  loader  is  orientated  at  0  degree  azimuth  and  0  degree 
elevation  with  reference  to  the  radar  direction  as  depicted  in  Figure  20: 


Figure  20  Orientation  of  CAD-Modeled  Backhoe  Loader  at  Aspect  Viewing 
Angles  of  Azimuth,  (0azimuth)=O°  and  Elevation,  (cpeievation)=0° 

The  radial  range  profiles  of  the  backhoe  loader  are  computed  using  MATLAB’s 
Fast  Fourier  Transform  of  the  associated  radar  backscattered  signals,  while  the  3-D  and  2 
-D  bi-range  profiles  are  computed  using  MATLAB’s  two-dimensional  (2-D)  discrete 
Fourier  Transform  of  the  third-order  cumulant  function  of  the  radar  backscatter  signals. 

Both  radial  range  and  bi-range  profiles  are  processed  over  a  frequency  bandwidth 
of  5.9058  GHz  (from  7.0472  to  12.953  GHz)  with  VV,  HH  and  HV  polarizations  for  the 
target’s  aspect  viewing  angle  of  azimuth,  (0azimuth)  =  0°  and  elevation,  (cpeievation)  =  0°. 
The  related  MATLAB  program,  “test  casel.m”  and  source  code  used  to  generate  this 
plot  can  be  found  in  Appendix  C. 
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Magnitude  of  Impulse  Response  Magnitude  of  Impulse  Response 


a.  Radial  Range  Profile  ( Gazimuth=0  <peievation=0  °) 

The  radial  range  profiles  of  the  backhoe  loader  with  VV,  HH  and  HV 
polarizations  are  illustrated  in  Figure  21,  22  and  23  respectively. 


Figure  2 1  Radial  Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 

0azimuth=O°  and  (peievation=0°  with  VV  Polarization 


Figure  22  Radial  Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 
0azimuth=O°  and  cpeievation=0°  with  HH  Polarization 
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Figure  23  Radial  Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 

0azimuth=O°  and  (peievation=0°  with  HV  Polarization 


b.  3-D  Bi-Range  Profile  ( Gazimuth=0  (pelevation  =0°) 

The  corresponding  3-D  bi-range  profiles  of  the  backhoe  loader  with  W, 
HH  and  HV  polarizations  are  shown  in  Figure  23,  24  and  25  respectively. 


Bi-range,  r2  (m)  Bi-range,  r,  (m) 

Figure  24  3-D  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 

0azimuth=O°  and  (peievation=0°  with  VV  Polarization 
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Magnitude  of  Bi-range  Response  Magnitude  of  Bi-range  Response 


Figure  25  3-D  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 

0azimuth=O°  and  (peievation=0°  with  HH  Polarization 


Bi-range,  r2  (m) 


Bi-range,  r,  (m) 


Figure  26  3-D  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 

0azimuth=O°  and  (peievation=0°  with  HV  Polarization 
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Bi-range,  r„  (m) 


c.  2-D  Characteristic  Bi-Range  Profile  ( 9azimuth=0  °,  <peievation=0  °) 

The  resulting  2-D  characteristic  bi-range  profiles  of  the  backhoe  loader 
with  VV,  HH  and  HV  polarizations  are  shown  in  Figure  25,  26  and  27  respectively. 


Figure  27  2-D  Characteristic  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect 

Angles  of  9azimuth=0o  and  (peievation=0o  with  VV  Polarization 
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Figure  28  2-D  Characteristic  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect 

Angles  of  9azimuth=0°  and  (peievation=0°  with  HH  Polarization 
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Bi-range,  r2  (m) 


Figure  29  2-D  Characteristic  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect 

Angles  of  9azimuth=0o  and  (peievation=0°  with  HV  Polarization 

3.  Test  Case  2:  Setup  and  Results  (0azimuth=9O°,  cpeievation=0°) 

In  Test  Case  2,  the  backhoe  loader  is  orientated  at  90  degree  azimuth  and  0  degree 
elevation  with  reference  to  the  radar  direction  as  depicted  in  Figure  30: 


Figure  30  Orientation  of  CAD-Modeled  Backhoe  Loader  at  Aspect  Viewing 
Angles  of  Azimuth,  (0azimuth)=90o  and  Elevation,  ((peievation)=0° 
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The  radial  range  profiles  of  the  backhoe  loader  are  computed  using  MATLAB’s 
Fast  Fourier  Transform  of  the  associated  radar  backscattered  signals,  while  the  3-D  and  2 
-D  bi-range  profiles  are  computed  using  MATLAB’s  two-dimensional  (2-D)  discrete 
Fourier  Transform  of  the  third-order  cumulant  function  of  the  radar  backscatter  signals. 

Both  radial  range  and  bi-range  profiles  are  processed  over  a  frequency  bandwidth 
of  5.9058  GHz  (from  7.0472  to  12.953  GHz)  with  VV,  HH  and  HV  polarizations  for  the 
target’s  aspect  viewing  angle  of  azimuth,  (©azimuth)  =  90°  and  elevation,  ((peievation)  =  0°. 
The  related  MATLAB  program,  “test_case2.m”  and  source  code  used  to  generate  this 
plot  can  be  found  in  Appendix  D. 


a.  Radial  Range  Profile  ( Ga-jmuth=90  °,  <peievation=0  °) 

The  radial  range  profiles  of  the  backhoe  loader  with  VV,  HH  and  HV 
polarizations  are  illustrated  in  Figure  31,  32  and  33  respectively. 


Figure  3 1  Radial  Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 
0azimuth=9O°  and  cpeievation=0o  with  VV  Polarization 
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Magnitude  of  Impulse  Response  Magnitude  of  Impulse  Response 


Figure  32  Radial  Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 
9azimuth=90°  and  (peievation=0°  with  HH  Polarization 


!  i  !  r 

i  i  i 

. . . .1 . ! . | ) 

.....1 . i . ! . 

_ i _ _ lJC 

_ 1 _ i _ 

Q _ I _ I  I  A  -'v _ *A. _ I _ v.  I _ I _ I _ 

-8  -6  -4  -2  0  2  4  6  8 

Range  (m) 


Figure  33  Radial  Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 
9azim.ith=90°  and  (peievation=9°  with  HV  Polarization 
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Magnitude  of  Bi-range  Response 


b.  3-D  Bi-Range  Profile  ( 0azimuth=9O  °,  elevation  =0°) 

The  corresponding  3-D  bi-range  profiles  of  the  backhoe  loader  with  W, 
HH  and  HV  polarizations  are  shown  in  Figure  34,  35  and  36  respectively. 


0  0 


Bi-range,  r-  (m) 


Bi-range,  ^  (m) 


Figure  34  3-D  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 

9azimuth=90°  and  cpeievation=0°  with  VV  Polarization 


0  0 


Bi-range,  r2  (m) 


Bi-range,  r,  (m) 


Figure  35  3-D  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 

9azimuth=90°  and  cpeievation=0o  with  HH  Polarization 
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Figure  36  3-D  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 

0azimuth=9O°  and  (peievation=0°  with  HV  Polarization 


c.  2-D  Characteristic  Bi-Range  Profile  ( 0azimuth=90  °,  (peievation=0  °) 
The  resulting  2-D  characteristic  bi-range  profiles  of  the  backhoe  loader 
with  VV,  HH  and  HV  polarizations  are  shown  in  Figure  37,  38  and  39  respectively. 
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Figure  37  2-D  Characteristic  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect 

Angles  of  0azimuth=9O°  and  cpeievation=0°  with  VV  Polarization 
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Figure  38  2-D  Characteristic  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect 

Angles  of  9azimuth=90o  and  cpeievation=0°  with  fffl  Polarization 


Figure  39  2-D  Characteristic  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect 

Angles  of  0azimuth=9O°  and  (peievation=0o  with  HV  Polarization 
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4.  Test  Case  3:  Setup  and  Results  (0azimuth=6O°,  (peievation=30°) 

In  Test  Case  3,  the  hackhoe  loader  is  orientated  at  60  degree  azimuth  and  30 
degree  elevation  with  reference  to  the  radar  direction  as  depicted  in  Figure  40: 


Figure  40  Orientation  of  CAD-Modeled  Backhoe  Loader  at  Aspect  Viewing 
Angles  of  Azimuth,  (0azimuth)=6O0  and  Elevation,  ((peievation)=30° 

The  radial  range  profiles  of  the  hackhoe  loader  are  computed  using  MATLAB’s 
Fast  Fourier  Transform  of  the  associated  radar  backscattered  signals,  while  the  3-D  and  2 
-D  bi-range  profiles  are  computed  using  MATLAB’s  two-dimensional  (2-D)  discrete 
Fourier  Transform  of  the  third-order  cumulant  function  of  the  radar  backscatter  signals. 

Both  radial  range  and  bi-range  profiles  are  processed  over  a  frequency  bandwidth 
of  5.9058  GHz  (from  7.0472  to  12.953  GHz)  with  VV,  HH  and  HV  polarizations  for  the 
target’s  aspect  viewing  angle  of  azimuth,  (0azimuth)  =  60°  and  elevation,  ((peievation)  =  30°. 
The  related  MATLAB  program,  “test_case3.m”  and  source  code  used  to  generate  this 
plot  can  be  found  in  Appendix  E. 
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Magnitude  of  Impulse  Response  Magnitude  of  Impulse  Response 


a.  Radial  Range  Profile  ( da-jmuth=60  °,  < peievation=30  °) 

The  radial  range  profiles  of  the  backhoe  loader  with  VV,  HH  and  HV 
polarizations  are  illustrated  in  Figure  41,  42  and  43  respectively. 


Figure  4 1  Radial  Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 
0azimuth=6O°  and  <peicvation:=30o  with  VV  Polarization 


Figure  42  Radial  Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 
0azimuth=6O°  and  (peievation=30°  with  HH  Polarization 
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Figure  43  Radial  Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 
0azimuth=6O°  and  (peievation=30°  with  HV  Polarization 


b.  3-D  Bi-Range  Profile  ( 0azimuth=60  (Pelevation  30  °) 

The  corresponding  3-D  bi-range  profiles  of  the  backhoe  loader  with  W, 
HH  and  HV  polarizations  are  shown  in  Figure  44,  45  and  46  respectively. 
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Figure  44  3-D  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 

0azimuth=6O°  and  cpeievation=30°  with  VV  Polarization 
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Figure  45  3-D  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 

0azimuth=6O°  and  cpeievation^O0  with  HH  Polarization 


Bi- range,  r2  (m)  Bi-range.  r,  (m) 

Figure  46  3-D  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect  Angles  of 

0azimuth=6O°  and  (peievation=30°  with  HV  Polarization 
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c.  2-D  Characteristic  Bi-Range  Profile  ( 9azimuth=90  °,  <peievation=0  °) 

The  resulting  2-D  characteristic  bi-range  profiles  of  the  backhoe  loader 
with  VV,  HH  and  HV  polarizations  are  shown  in  Figure  47,  48  and  49  respectively. 


Figure  47  2-D  Characteristic  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect 

Angles  of  9azimuth=60o  and  cpeievation^O0  with  VV  Polarization 


Figure  48  2-D  Characteristic  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect 

Angles  of  9azimuth=60°  and  cpeievation^O0  with  HH  Polarization 
52 


Figure  49  2-D  Characteristic  Bi-Range  Profile  of  Backhoe  Loader  at  Aspect 

Angles  of  0azimuth=6O°  and  (peievation=30°  with  HV  Polarization 

B.  ANALYSIS  OF  “BACKHOE  DATA  DOME”  EXPERIMENTATION 

In  studying  and  comparing  the  various  radial  range  and  bi-range  profiles 
computed  using  the  “backhoe  data  dome”  experimentation  for  Test  Cases  1,  2  and  3, 
observations  and  analysis  relating  to  the  uses  of  radial  range  profiles,  polarization  effects 
and  bi-range  profiles  for  target  recognition  can  be  rationalized  and  are  discussed  in  the 
following  sections. 

1.  Radial  Range  Profiles 

The  magnitude  of  the  impulse  response  versus  the  range  plots  of  the 
experimentation  results  provide  a  good  representation  of  the  one-dimensional  radial 
profile  of  the  backhoe  loader  from  different  aspect  angles.  It  can  be  observed  that 
locations  of  the  peak  responses  appear  to  coincide  and  correspond  to  the  relative 
positions  of  dominant  individual  scattering  centers  or  subcomponents,  such  as  the 
operator  cab,  engine  compartment,  and  loader  bucket,  etc.,  that  present  a  large  RCS  and 
consequently  a  strong  backscattered  signal  to  the  radar  system  as  illustrated  in  Figure  50. 
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Figure  50  Comparison  of  Radar  Range  Profile  with  Actual  Backhoe  Loader 

While  it  may  seem  that  the  radial  range  profile  can  possibly  provide  a  measure  of 
the  length  dimension  of  a  target  of  interest  along  its  line-of-sight  from  the  radar  system, 
this  inference  is  usually  erroneous  and  not  equivalent  to  the  actual  physical  target’s 
length  due  to  two  effects,  as  mentioned  in  Tait  (2005).  Instead  of  directly  reflecting  the 
radar  signal  back  toward  the  radar  system,  multiple  scatterers  comprising  the  target  and 
their  interactions  can  result  in  indirect  reflection  and  delay  the  round-trip  propagation 
time,  which  then  produces  a  radial  range  profile  that  is  longer  than  the  target. 
Conversely,  the  other  effect  results  in  a  shorter  radial  range  profile  compared  to  the  actual 
physical  target’s  length.  This  disparity  is  due  to  the  shadowing  or  masking  of  one  part 
the  target  by  another  part.  When  a  certain  part  of  a  target  is  masked  at  a  particular  aspect 
angle,  the  obscuration  reduces  the  electromagnetic  energy  that  scattered  off  of  it,  and 
significantly  reduces  its  contribution  to  the  radial  range  profile. 
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For  example,  as  also  observed  in  Figure  50,  the  dipperstick  of  the  backhoe  loader 
is  obscured  by  the  engine  compartment  and  operator  cab,  which  then  reduces  its 
corresponding  response  in  the  radial  range  profile  and  results  in  an  apparently  shorter 
backhoe  loader. 

2.  Polarization  Effects 

The  polarization  of  a  radar  signal  is  defined  as  the  direction  of  the  electric  E-field 
as  the  electromagnetic  wave  propagates.  When  a  target  of  interest  is  illuminated  by  a 
linearly  polarized  signal,  two  echo  signals  are  reflected  back  to  the  radar  system;  one 
with  the  same  polarization  as  the  transmitted  signal  (co-polarized  signal)  and  the  other 
with  an  orthogonal  polarization  (cross-polarized  signal).  The  basis  of  using  polarization 
effects  for  target  recognition  is  to  extract  the  differences  in  the  radar  backscattered 
signals  from  a  target  that  is  illuminated  with  different  polarization.  Although  this  method 
seems  to  be  promising  for  discriminating  simple  geometrically  shaped  targets,  such  as 
spheres,  comer  reflectors  and  straight  lengths  of  wire,  etc,  it  is  not  as  forthcoming  for 
practical  complex  targets,  such  as  aircraft  and  vehicles. 

For  example,  it  is  noted  that  the  3-D  and  2-D  bi-range  profiles  of  the  backhoe 
loader  for  the  different  test  cases  are  not  affected  by  the  polarization  of  the  radar  signal  as 
shown  by  the  similarities  between  the  different  plots  with  VV  and  HH  polarizations. 

This  is  probably  due  to  the  multiple  scatterers  or  subcomponents  interaction 
within  the  radar  resolution  cell  that  collectively  resulted  in  a  composite  radar 
backscattered  signal  with  a  different  and  independent  polarization  as  the  transmitted 
signal.  In  addition,  it  is  also  observed  from  the  3-D  and  2-D  bi-range  profiles  with  HV 
polarization  that  the  peak  responses  from  multiple  scatterers’  interactions  seems  to  be 
accentuated  by  the  cross  polarization  effect. 

3.  Characteristic  Bi-Range  Profile 

Most  importantly,  it  is  obvious  from  the  experimentation  results  of  the  different 
test  cases  that  there  exists  significant  peak  responses  in  the  3-D  and  2-D  bi-range  profiles, 
implying  strong  additive  multiple  interactions  between  the  scatterers  or  subcomponents 
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of  the  backhoe  loader,  that  are  characteristic  of  their  bi-range  shapes  and  locations.  The 
positions  of  the  multiple  scattering  interaction  terms,  relative  to  the  direct  specular 
scattering  tenns,  correspond  directly  to  the  physical  separation  distances  between  the 
major  individual  scatterers  of  the  backhoe  loader. 

Therefore,  these  characteristic  bi-range  profiles  can  be  drawn  on,  to  highlight  and 
map  the  unique  multiple  scattering  interaction  features  and  natures  of  the  target  which 
may  then  be  used  to  compare  and  match  against  a  reference  database  that  contains  bi¬ 
range  signatures  of  different  target  types  of  potential  interest  and  at  all  aspect  angles  of 
interest  for  target  recognition  and  identification. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


Among  the  detection  systems  under  development  in  the  DD(X)  program  is 
the  Multi-Function  Radar  (MFR).  The  luturc  MFR  will  provide  Non- 
Cooperative  Target  Recognition  (NCTR)  and  Kill  Assessment. 

—  Rear  Admiral  Philip  Balisle,  U.S.  Navy  (2002) 


A.  SUMMING  UP  OF  THESIS  RESEARCH 

Radar-based  NCTR  will  be  a  key  capability  and  cardinal  requirement  in  future 
radar  systems  equipped  onboard  all  high-valued  military  combat  platforms  and  assets, 
such  as  fighter  aircraft  and  destroyers,  etc.  While  there  are  many  conventional 
techniques  for  the  incorporation  of  NCTR  into  the  next-generation  radar  systems,  this 
thesis  endeavors  to  expand  on  the  knowledge  and  utility  of  the  bi-spectral  method  for 
effective  and  robust  target  recognition  and  identification. 

The  first  objective  of  this  thesis  was  to  study  and  mathematically  model  the  direct 
and  indirect  multiple  interactions  of  scattering  mechanisms  of  radar  backscattered  signal. 
The  fundamentals  of  radar  basic  operation,  electromagnetic  wave  propagation  and 
scattering  mechanisms  of  electromagnetic  waves,  when  a  radar  signal  intersects  any 
target  in  its  propagation  path,  were  first  reviewed.  This  then  provided  an  insight  into  the 
relationships  between  the  target’s  radar  signature  and  its  physical  size,  geometry,  shape, 
orientation  and  material  composition,  which  were  the  basis  in  formulating  the 
mathematical  model  of  the  radar  backscattered  signal. 

A  mathematical  model  which  represents  the  direct  specular  and  indirect  multiple 
scattering  interactions  of  the  backscattered  radar  signal  was  developed  and  subsequently 
verified  using  an  example  of  a  simple  three-point  scatterers  target  and  the  MATLAB 
program  “backscatter  signal.m”  to  plot  the  radial  range  profile. 

The  second  objective  of  this  thesis  was  to  develop  a  radar  NCTR  technique  of 
using  the  bi-spectral  signatures  of  backscattered  radar  signals  for  characterization  of  a 
target’s  multiple  interactions  features.  The  definitions  and  mathematical  derivations  of 
the  bi-spectrum,  bi-time  and  bi-range  of  the  radar  backscattered  signal  were  introduced, 
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together  with  the  concept  of  using  the  bi-spectral  signatures,  also  known  as  the 
characteristic  bi-range  profiles  of  backscattered  radar  signal  as  a  complementary  add-on 
to  augment  the  radial  range  profile  for  target  recognition  and  identification. 

The  characteristic  bi-range  profile  was  then  studied  and  verified  through  the  same 
example  of  a  simple  3-points  scatterer  target  and  the  MATLAB  program,  “bispectral.m” 
was  developed  to  compute  and  plot  the  bi-range  profiles  of  the  third-order  cumulant  or 
triple  correlation  of  the  frequency-domain  radar  backscattered  signal. 

To  further  investigate  the  potential  of  using  bi-spectra  technique  for  radar  target 
recognition,  experimentation  based  on  the  AFRL,  “Backhoe  Data  Dome”  was  conducted. 
The  results  and  findings  for  three  different  test  cases  were  analyzed  and  found  to  be 
promising. 

B.  RECOMMENDATIONS  FOR  FUTURE  WORK 

The  novel  idea  of  using  the  bi-spectral  method  for  radar  target  recognition,  based 
on  the  characterization  of  a  target’s  multiple  scatterer  interactions  tenn,  is  a  relatively 
recent  development,  considering  that  the  concept  was  initially  conceived  and  proposed  as 
a  basis  of  recognition  by  Jounty  &  Walton  (1990)  during  the  early  nineties.  As  such, 
many  aspects  of  this  proposed  bi-spectal  NCTR  technique  can  be  further  researched, 
developed  and  verified.  The  following  sections  recommend  two  areas  for  which  future 
research  and  investigation  can  be  considered. 

1.  Experimentation  Based  on  Actual  Radar  Signature  Measurements  of 
Real  Complex  Target 

For  the  current  thesis,  due  to  the  restricted  access  to  classified  infonnation,  the 
experimentation  to  examine  the  potential  of  using  the  characteristic  bi-range  profile  for 
target  recognition  and  identification  is  based  on  the  AFRL’s  public  release  of  the 
“Backhoe  Data  Dome,  Version  1.0,”  which  consists  of  simulated  wideband  (7-13  GHz), 
full  polarization,  complex  backscatter  data  from  a  backhoe  loader  generated  through  the 
use  of  a  CAD  model  in  free  space  (Lin  &  Naidu,  2004).  While  this  set  of  data  facilitated 
the  verification  of  the  concept  of  using  a  bi-spectral  method  for  target  recognition,  the 
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experimentation  lacks  fidelity  in  the  consideration  of  the  radar  bi-spectral  signatures  and 
scattering  mechanisms  from  a  real  complex  target  in  an  environment  surrounded  with 
clutter  and  noise. 

Given  proper  security  clearance,  more  comprehensive  experimentation,  based  on 
classified  data  sets  of  actual  radar  backscattered  signal  measurements  from  real  complex 
targets  can  be  examined  and  analyzed  in  future  research  to  better  understand  the 
robustness  of  the  bi-spectral  method  and  its  clutter  and  noise  suppression  performance. 

2.  Bi-spectral  Waveform  Design 

Another  area  of  research  interest  that  can  be  further  investigated  and  verified  for 
future  work  is  the  topic  of  bi-spectral  wavefonn  design.  The  design  of  the  radar 
waveform  is  crucial  in  the  measurement  of  the  radar  backscattered  signals  required  to 
generate  the  bi-spectral  profile.  In  order  to  generate  bi-range  profiles,  which  are 
compatible  with  the  reference  signature  database  of  different  target  types  of  potential 
interest  for  target  recognition,  the  radar  waveform  needs  to  provide  high  range  resolution 
measurements  and  extract  discriminating  multiple  scatterers  interaction  features  of  the 
target. 

While  the  experimentation  results  from  the  current  thesis  seem  to  suggest  that  the 
use  of  a  cross  polarized  signal  can  possibly  accentuate  the  peak  bi-range  profile  responses 
from  multiple  scatterers’  interactions,  more  in-depth  investigation  and  studies  have  to  be 
performed  in  order  to  understand  and  validate  the  polarization  characteristics  of  the 
multiple  scatterers’  interactions  terms.  Apart  from  the  polarization  of  the  radar  signal, 
other  radar  waveform  design  parameters  like  the  transmitted  frequency,  pulse  width, 
pulse  repetition  frequency,  frequency  bandwidth  and  modulation,  etc.,  can  also  be 
examined  in  future  research  work  in  order  to  establish  an  optimal  bi-spectral  waveform. 
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APPENDIX 


MATLAB  SOURCE  CODES 


A.  SOURCE  CODE  FOR  RADAR  BACKSCATTERED  SIGNAL 

% - 

%  backscatter_signal.m 

% 

%  This  program  computes  the  radar  backscattered  signal  and 
%  generates  the  radial  range  profile  of  n  points  scatterers 

% 

%  Written  by  Y eo,  Jiunn  W ah 

% - 

clear  all 


n  =  3; 

% 

Number  of  point  scatter 

f=  5:0.01:6; 

% 

Frequency  band  in  GHz 

A  =  [1  1  1]; 

% 

Amplitude  of  the  scattered  signal 

R  =  [0  1.5  3.6]; 

% 

Range  between  the  Radar  and  point  scatter 

s  =  0; 

% 

Initialization  of  counter 

c  =  3e8; 

% 

Speed  of  light 

bandwidth  =  f(length(f))-f(l); 

%  Direct  specular  scattering  term 
for  count  1=1  :n 

s=s+A(countl).*exp(i.*2.*pi.*2*R(countl).*le9.*f./c);%-rand*2*pi); 

end 


%  Indirect  multiple  scattering  interaction  term 
for  count2=l  :n 
for  count3=l:n 
if  count3~=count2 

s=s+A(count2)*A(count3).*exp(i.*2.*pi.*2*(R(count2)+R(count3)).*le9.*f./c); 

end 

end 

end 


%  Plotting  of  Radial  Range  Profile 

z=fft(s)./length(s); 

figure 

plot(linspace(-length(z)/2,length(z)/2,length(z))*c/2/bandwidth/le9,  fftshift(abs(z))); 
xlabel('Range  (m)') 
ylabel('Magnitude  Impulse  Response') 
grid  on 
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B. 


SOURCE  CODE  FOR  GENERATING  THE  THIRD-ORDER  CUMULANT 
AND  BI-SPECTRAL  PLOTS 


% - 

%  bispectral.m 

% 

%  This  program  generates  the  3-D  and  2-D  bi-range  profile  of 

%  a  radar  backscattered  signal 

% 

%  Written  by  Yeo,  Jiunn  Wah 

% - 

clear  all 


%  Call  program  to  compute  the  radar  backscattered  signal  and 
%  generates  the  radial  range  profile  of  n  points  scatterer 
backscattersignal; 

%  Function  to  generate  the  third  order  cumulant  of  the  radar  backscattered  signal 
[ex]  =  cum3(s,(length(s)-l)); 

%  Perform  Fourier  Transform  of  the  third  order  cumulant  to  generate  the  bispectrum 
bisp  =  fftshift(fft2(ifftshift(cx))); 

%  Setting  the  plot  limits 

x  =  linspace(-length(bisp)/4,length(bisp)/4,length(bisp))*c/2/bandwidth/le9; 
y  =  linspace(-length(bisp)/4,length(bisp)/4,length(bisp))*c/2/bandwidth/le9; 
ref  =  find(x==0); 

%  Plotting  of  the  3-D  bi-range  profile 
figure 

mesh(x(ref:length(x)),y(ref:length(y)),abs(bisp(ref:length(bisp),(ref  length(bisp))))); 
xlabel('Bi-range,  r_l  (m)') 
ylabel('Bi-range,  r_2  (m)') 
zlabel('Magnitude  of  Bi-range  Response') 

%  Plotting  of  the  2-D  Characteristic  bi-range  profile 
figure 

contourf(x(ref:length(x)),y(ref:length(y)),abs(bisp(ref:length(bisp),(ref:length(bisp)))),  [2e4  1.2e5  1.4e5 

1.6e5  1.8e5  2e5|): 

xlabel('Bi-range,  r_l  (m)') 

ylabel('Bi-range,  r_2  (m)') 

grid  on 

colorbar 


% - 

%  cum3.m 

% 

%  This  function  generates  the  third  order  cumulant  of  the  radar  backscattered  signal 

% 

%  Written  by  Y eo,  Jiunn  W ah 

% 

%  ex  :  the  cross  cumulant  matrix  C(m,n)  of  the  radar  backscattered  signal 

%  signal  :  the  radar  backscattered  signal  to  be  processed 

% 

%  signal  length  :  the  maximum  argument  of  the  cross  cumulant 

% - 
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function  [cx]  =  cum3(signal,signaMength) 

%  Assignment  and  pre-conditioning  of  signal 

signal_l=conj(signal(:)); 

signal_2=signal(:); 

signal_3=signal(:); 

N  =  length(signal_l); 

%  Triple  correlation  of  signal 

XI  =  repmat(signal_l,l,2*signaMength+l); 

C2=[zeros(signal_length,l);signal_2(l:signal_length+l)]; 
R2=[signal_2(signal_length+1:N).'  zeros(l,signaMength)]; 
X2=hankel(C2,R2); 

C3=[zeros(signal_length,l);signal_3(l:signal_length+l)]; 
R3=[signal_3(signal_length+1:N).'  zeros(l,signaMength)]; 
X3=hankel(C3  ,R3 

cx=X2*  (X 1 .  *X3  )/N ; 


C.  SOURCE  CODE  FOR  TEST  CASE  1 

% - 

%  test_casel.m 

% 

%  This  program  computes  the  third  order  cumulant  and  bispectrum  of  the  “backhoe 
%  data  dome”  radar  backscattered  signals  for  a  target  orientation  of  azimuth  angle, 

%  (^azimuth)  =  0°  and  elevation  anlge,  ((peievation)  =  0°. 

% 

%  Both  radial  range  and  bi-range  profiles  are  processed  and  plotted  over  a  frequency 
%  bandwidth  of  5.9058  GHz  (from  7.0472  to  12.953  GHz)  with  W,  HH  and  HV 
%  polarizations. 

% 

%  Written  by  Yeo,  Jiunn  Wah 

% - 

clear  all; 

%  Load  “Backhoe  Data  Dome”  data 

load('C:\DOCUME~l\J&C\MYDOCU~l\NPS\Thesis\BACKHO~l\2D_CHA~l\2D_K_S~l\backhoe_elOO 

0_az350tol00.mat') 


c  =  3e8; 

% 

Speed  of  light 

bandwidth  =  5.9058; 

% 

Stimulated  frequency  bandwidth 

sVV  0  =  VV(:,141); 
zVV_0  =  fft(sVV_0); 

% 

VV  polarized  radar  signal  at  0azimuth  =  0°  and 

^Pelevation  0 

sHH  0  =  HH(:,141); 
zHHO  =  fft(sHH_0); 

% 

HH  polarized  radar  signal  at  0azimuth  =  0°  and 

^Pelevation  —  0 

sHV_0  =  HV(:,141); 

% 

HV  polarized  radar  signal  at  0azimuth  =  0°  and 

^Pelevation  —  0 
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zHVO  =  fft(sHV_0); 


%  Plotting  of  Radial  Range  Profile  with  VV  polarization 
figure 

plot((-length(zVV_0)/2:length(zVV_0)/2-l)*c/2/bandwidth/le9,  (abs(fftshift(zVV_0)))) 
xlabel('Range  (m)') 

ylabel('Magnitude  of  Impulse  Response') 

title('Radial  Range  Profile  for  \theta_a_z  =  0  deg  and  \theta_e_l  =  0  deg  with  VV  Polarization') 
grid  on 


%  Plotting  of  Radial  Range  Profile  with  HH  polarization 
figure 

plot((-length(zHH_0)/2:length(zHH_0)/2-l)*c/2/bandwidth/le9,  fftshift(abs((zHH_0)))) 
xlabel('Range  (m)') 

ylabel('Magnitude  of  Impulse  Response') 

title('Radial  Range  Profile  for  \theta_a_z  =  0  deg  and  \theta_e_l  =  0  deg  with  HH  Polarization') 
grid  on 


%  Plotting  of  Radial  Range  Profile  with  HV  polarization 

figure 

plot((-length(zHV_0)/2:length(zHV_0)/2-l)*c/2/bandwidth/le9,  (abs(fftshift(zHV_0)))) 
xlabel('Range  (m)') 

ylabel('Magnitude  of  Impulse  Response') 

title('Radial  Range  Profile  for  \theta_a_z  =  0  deg  and  \theta_e_l  =  0  deg  with  HV  Polarization') 
grid  on 


%  Computation  of  third  order  cumulant  and  bispectrum 
[cxVVO]  =  cum3  ( s  V  V_0,  ( length)  s  V  V_0)- 1 )) ; 
bispVVO  =  fftshift(fft2(ifftshift(cxVV_0))); 

[cxHHO]  =  cum3(sHH_0, (length)  sHHO)- 1 )); 
bispHH_0  =  fftshift(fft2(ifftshift(cxHH_0))); 

[cxHVO]  =  cum3(sHV_0, (length)  sHVO)- 1 )); 
bispHV_0  =  fftshift(fft2(ifftshift(cxHV_0))); 

%  Plotting  of  the  3-D  bi-range  profile  with  VV  polarization 
figure 

xVV_0=linspace(-length(bispVV_0)/4,length(bispVV_0)/4,length(bispVV_0))*c/2/bandwidth/le9; 
yVV_0  =  linspace(-length(bispVV_0)/4,length(bispVV_0)/4,length(bispVV_0))*c/2^andwidth/le9; 
refVVO  =  fmd(xVV_0==0); 
re  fend  VV_0  =  find(xVV_0>=4.5  &  xVV_0<=4.52); 

mesh(xVV_0(refVV_0:refendVV_0),yVV_0(refVV_0:refendVV_0),abs(bispVV_0(refVV_0:refendVV_0, 

(re  fV  V_0 :  refend  V  V_0 ) ) ) ) ; 

xlabel('Bi-range,  r_l  (m)') 

ylabel('Bi-range,  r_2  (m)') 

zlabel('Magnitude  of  Bi-range  Response') 

title('3-D  Bi-range  Profile  for  \theta_a_z  =  0  deg  and  \theta_e_l  =  0  deg  with  VV  Polarization') 

%  Plotting  of  the  3-D  bi-range  profile  with  HH  polarization 
figure 

xHH_0  =  linspace(-length(bispHH_0)/4,length(bispHH_0)/4,length(bispHH_0))*c/2/bandwidth/le9; 
yHHO  =  linspacc(-lcngth(bispHH_0)4.length(bispHH_0)4.length(bispHH_0))*c/2/bandwidthl  c9; 
refHHO  =  find(xHH_0==0); 
refendHH  O  =  find(xHH_0>=4.5  &  xHH_0<=4.52); 
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mesh(xHH_0(refHH_0:refendHH_0),yHH_0(refHH_0:refendHH_0),abs(bispHH_0(refHH_0:refendHH_0, 

(refHH_0:refendHH_0)))); 

xlabel('Bi-range,  r_l  (m)') 

ylabel('Bi-range,  r_2  (m)') 

zlabel('Magnitude  of  Bi-range  Response') 

title('3-D  Bi-range  Profile  for  \theta_a_z  =  0  deg  and  \theta_e_l  =  0  deg  with  HH  Polarization') 

%  Plotting  of  the  3-D  bi-range  profile  with  HV  polarization 
figure 

xHVO  =  linspace(-length(bispHV_0)/4,length(bispHV_0)/4,length(bispHV_0))*c/2/bandwidth/le9; 
yHV_0  =  hnspace(-length(bispHV_0)/4,length(bispHV_0)/4,length(bispHV_0))*c/2/bandwidth/le9; 
refHVO  =  fmd(xHV_0>=-2.51  &  xHV_0<=-2.5); 
refendHVO  =  find(xHV_0>=2.5  &  xHV_0<=2.51); 

mesh(xHV_0(refHV_0:refendHV_0),yHV_0(refHV_0:refendHV_0),abs(bispHV_0(refHV_0:refendHV_0, 

(refHV_0:refendHV_0)))); 

xlabel('Bi-range,  r_l  (m)') 

ylabel('Bi-range,  r_2  (m)') 

zlabel('Magnitude  of  Bi-range  Response') 

title('3-D  Bi-range  Profile  for  \theta_a_z  =  0  deg  and  \theta_e_l  =  0  deg  with  HV  Polarization') 

%  Plotting  of  the  2-D  characteristic  bi-range  profile  with  VV  polarization 
figure 

contour(xVV_0(refVV_0:refendVV_0),yVV_0(refVV_0:refendVV_0),abs(bispVV_0(refVV_0:refendVV 
_0 ,  (re  fV  V_0 :  re  fend  V  V_0 )) ) ,  [2e7  4e7  6e7  8e7  le8  1.2e8  1.4e8  1.6e8  1.8e8  2e8]); 
xlabel('Bi-range,  r_l  (m)') 
ylabel('Bi-range,  r_2  (m)') 

title('2-D  Characteristic  Bi-range  Profile  for  \theta_a_z  =  0  deg  and  \theta_e_l  =  0  deg  with  VV 

Polarization') 

grid  on 

colorbar 


%  Plotting  of  the  2-D  characteristic  bi-range  profile  with  HH  polarization 
figure 

contour(xHH_0(refHH_0:refendHH_0),yHH_0(refHH_0:refendHH_0),abs(bispHH_0(refHH_0:refendHH 
_0 , (re fHHO : re fendHHO )) ) ,  [2e7  4e7  6e7  8e7  le8  1.2e8  1.4e8  1.6e8  1.8e8  2e8]); 
xlabel('Bi-range,  r_l  (m)') 
ylabel('Bi-range,  r_2  (m)') 

title('2-D  Characteristic  Bi-range  Profile  for  \theta_a_z  =  0  deg  and  \theta_e_l  =  0  deg  with  HH 

Polarization') 

grid  on 

colorbar 


%  Plotting  of  the  2-D  characteristic  bi-range  profile  with  HV  polarization 
figure 

contour(xHV_0(refHV_0:refendHV_0),yHV_0(refHV_0:refendHV_0),abs(bispHV_0(refHV_0:refendHV 
_0 , (re fH V  O : re fendH V  O )) ) ,  [3  5  10  15  20  25  30]); 
xlabel('Bi-range,  r_l  (m)') 
ylabel('Bi-range,  r_2  (m)') 

title('2-D  Characteristic  Bi-range  Profile  for  \theta_a_z  =  0  deg  and  \theta_e_l  =  0  deg  with  HV 

Polarization') 

grid  on 

colorbar 
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D.  SOURCE  CODE  FOR  TEST  CASE  2 

% - 

%  test_case2.m 

% 

%  This  program  computes  the  third  order  cumulant  and  bispectrum  of  the  “backhoe 
%  data  dome”  radar  backscattered  signals  for  a  target  orientation  of  azimuth  angle, 

%  (^azimuth)  =  90°  and  elevation  anlge,  (<peievation)  =  0°. 

% 

%  Both  radial  range  and  bi-range  profiles  are  processed  and  plotted  over  a  frequency 

%  bandwidth  of  5.9058  GHz  (from  7.0472  to  12.953  GHz)  with  VV,  HH  and  HV 

%  polarizations. 

% 

%  Written  by  Y eo,  Jiunn  W ah 

% - 

clear  all; 


%  Load  “Backhoe  Data  Dome”  data 

load('C:\DOCUME~l\J&C\MYDOCU~l\NPS\Thesis\BACKHO~l\2D_CHA~l\2D_K_S~l\backhoe_elOO 

0_az350tol00.mat') 


c  =  3e8; 

% 

Speed  of  light 

bandwidth  =  5.9058; 

% 

Stimulated  frequency  bandwidth 

sVV  90  =  VV(:,1401); 
zVV_90  =  fft(sVV_90); 

% 

VV  polarized  radar  signal  at  0azimuth  =  90°  and 

elevation  0 

sHH  90  =  HH(:,1401); 
zHH_90  =  fft(sHH_90); 

% 

HH  polarized  radar  signal  at  0azimuth  =  90°  and 

^Pelevation  0 

sHV  90  =  HV(:,1401); 
zHV_90  =  fft(sHV_90); 

% 

HV  polarized  radar  signal  at  0azimuth  =  90°  and 

^Pelevation  —  0 

%  Plotting  of  Radial  Range  Profile  with  VV  polarization 
figure 

plot((-length(zVV_90)/2:length(zVV_90)/2-l)*c/2/bandwidth/le9,  (abs(fftshift(zVV_90)))) 
xlabel('Range  (m)') 

ylabel('Magnitude  of  Impulse  Response') 

title('Radial  Range  Profile  for  \theta_a_z  =  90  deg  and  \theta_e_l  =  0  deg  with  VV  Polarization') 
grid  on 


%  Plotting  of  Radial  Range  Profile  with  HH  polarization 
figure 

plot((-length(zHH_90)/2:length(zHH_90)/2-l)*c/2/bandwidth/le9,  fftshift(abs((zHH_90)))) 
xlabel('Range  (m)') 

ylabel('Magnitude  of  Impulse  Response') 

title('Radial  Range  Profile  for  \theta_a_z  =  90  deg  and  \theta_e_l  =  0  deg  with  HH  Polarization') 
grid  on 


%  Plotting  of  Radial  Range  Profile  with  HV  polarization 
figure 

plot((-length(zHV_90)/2:length(zHV_90)/2-l)*c/2/bandwidth/le9,  (abs(fftshift(zHV_90)))) 
xlabel('Range  (m)') 

ylabel('Magnitude  of  Impulse  Response') 

title('Radial  Range  Profile  for  \theta_a_z  =  90  deg  and  \theta_e_J  =  0  deg  with  HV  Polarization') 
grid  on 
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%  Computation  of  third  order  cumulant  and  bispectrum 
[cxVV_90]  =  cum3(sW_90,(length(sW_90)-l)); 
bispVV_90  =  fftshift(fft2(ifftshift(cxVV_90))); 

[cxHH_90]  =  cum3(sHH_90,(length(sHH_90)-l)); 
bispHH_90  =  fftshift(fft2(ifftshift(cxHH_90))); 

[cxHV_90]  =  cum3(sHV_90,(length(sHV_90)-l)); 
bispHV_90  =  fftshift(fft2(ifftshift(cxHV_90))); 

%  Plotting  of  the  3-D  bi-range  profile  with  VV  polarization 
figure 

xVV_90  =  linspace(-length(bispVV_90)/4,length(bispVV_90)/4,length(bispVV_90))*c/2/bandwidth/le9; 
yVV_90  =  linspace(-length(bispVV_90)/4,length(bispVV_90)/4,length(bispVV_90))*c/2/bandwidth/le9; 
refVV_90  =  find(xW_90==0); 
refendVV_90  =  find(xVV_90>=2  &  xVV_90<=2.0085); 

mesh(xW_90(refW_90:refendW_90),yW_90(refW_90:refendW_90),abs(bispW_90(refW_90:ref 

endVV_90,(refVV_90:refendVV_90)))); 

xlabel('Bi-range,  r_l  (m)') 

ylabel('Bi-range,  r_2  (m)') 

zlabel('Magnitude  of  Bi-range  Response') 

title('3-D  Bi-range  Profile  for  \theta_a_z  =  90  deg  and  \theta_e_l  =  0  deg  with  VV  Polarization') 

%  Plotting  of  the  3-D  bi-range  profile  with  HH  polarization 
figure 

xHH_90  =  linspace(-length(bispHH_90)/4,length(bispHH_90)/4,length(bispHH_90))*c/2/bandwidth/le9; 
yHH_90  =  linspace(-length(bispHH_90)/4,length(bispHH_90)/4,length(bispHH_90))*c/2/bandwidth/le9; 
refHH_90  =  find(xHH_90==0); 
refendHH_90  =  find(xHH_90>=2  &  xHH_90<=2.0085); 

mesh(xHH_90(refHH_90:refendHH_90),yHH_90(refHH_90:refendHH_90),abs(bispHH_90(refHH_90:ref 

endHH_90,(refHH_90:refendHH_90)))); 

xlabel('Bi-range,  r_l  (m)') 

ylabel('Bi-range,  r_2  (m)') 

zlabel('Magnitude  of  Bi-range  Response') 

title('3-D  Bi-range  Profile  for  \theta_a_z  =  90  deg  and  \theta_e_l  =  0  deg  with  HH  Polarization') 

%  Plotting  of  the  3-D  bi-range  profile  with  HV  polarization 
figure 

xHV_90  =  linspace(-length(bispHV_90)/4,length(bispHV_90)/4,length(bispHV_90))*c/2/bandwidth/le9; 
yHV_90  =  linspace(-length(bispHV_90)/4,length(bispHV_90)/4,length(bispHV_90))*c/2/bandwidth/le9; 
refHV_90  =  find(xHV_90>=-2.01  &  xHV_90<=-2.0); 
refendHV_90  =  find(xHV_90>=2.0  &  xHV_90<=2.01); 

mesh(xHV_90(refHV_90:refendHV_90),yHV_90(refHV_90:refendHV_90),abs(bispHV_90(refHV_90:ref 

endHV_90,(refHV_90:refendHV_90)))); 

xlabel('Bi-range,  r_l  (m)') 

ylabel('Bi-range,  r_2  (m)') 

zlabel('Magnitude  of  Bi-range  Response') 

title('3-D  Bi-range  Profile  for  \theta_a_z  =  90  deg  and  \theta_e_l  =  0  deg  with  HV  Polarization') 

%  Plotting  of  the  2-D  characteristic  bi-range  profile  with  VV  polarization 
figure 

contour(xVV_90(refVV_90:refendVV_90),yVV_90(refVV_90:refendVV_90),abs(bispVV_90(refVV_90:r 
efendVV_90,(refVV_90:refendVV_90))),  [2.5e9  5e9  7.5e9  lelO  1.5el0  2el0  2.5el0]); 
xlabel('Bi-range,  r_l  (m)') 
ylabel('Bi-range,  r_2  (m)') 
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title('2-D  Characteristic  Bi-range  Profile  for  \theta_a_z  =  90  deg  and  \theta_e_l  =  0  deg  with  VV 

Polarization') 

grid  on 

colorbar 


%  Plotting  of  the  2-D  characteristic  bi-range  profile  with  HH  polarization 
figure 

contour(xHH_90(refHH_90:refendHH_90),yHH_90(refHH_90:refendHH_90),abs(bispHH_90(refHH_90:r 
efendHH_90,(refHH_90:refendHH_90))),  [2.5e9  5e9  7.5e9  lelO  1.5el0  2el0  2.5el0]); 
xlabel('Bi-range,  r_l  (m)') 
ylabel('Bi-range,  r_2  (m)') 

title('2-D  Characteristic  Bi-range  Profile  for  \theta_a_z  =  90  deg  and  \theta_e_l  =  0  deg  with  HH 

Polarization') 

grid  on 

colorbar 


%  Plotting  of  the  2-D  characteristic  bi-range  profile  with  HV  polarization 
figure 

contour(xHV_90(refHV_90:refendHV_90),yHV_90(refHV_90:refendHV_90),abs(bispHV_90(refHV_90:r 
efendHV_90,(refHV_90:refendHV_90))),  [18  20  40  60  80  100  120  140  160  180]); 
xlabel('Bi-range,  r_l  (m)') 
ylabel('Bi-range,  r_2  (m)') 

title('2-D  Characteristic  Bi-range  Profile  for  \theta_a_z  =  90  deg  and  \theta_e_l  =  0  deg  with  HV 

Polarization') 

grid  on 

colorbar 

E.  SOURCE  CODE  FOR  TEST  CASE  3 

% - 

%  test_case3.m 

% 

%  This  program  computes  the  third  order  cumulant  and  bispectnim  of  the  “backhoe 
%  data  dome”  radar  backscattered  signals  for  a  target  orientation  of  azimuth  angle, 

%  (dazimuth)  =  60°  and  elevation  anlge,  ((pcievation)  =  30°. 

% 

%  Both  radial  range  and  bi-range  profiles  are  processed  and  plotted  over  a  frequency 

%  bandwidth  of  5.9058  GHz  (from  7.0472  to  12.953  GHz)  with  VV,  HH  and  HV 

%  polarizations. 

% 

%  Written  by  Y eo,  Jiunn  W ah 

% - 

clear  all; 


%  Load  “Backhoe  Data  Dome”  data 

load('C:\DOCUME~l\J&C\MYDOCU~l\NPS\Thesis\BACKHO~l\2D_CHA~l\2D_K_S~l\backhoe_el03 

0_az350tol00.mat') 


c  =  3e8; 

% 

Speed  of  light 

bandwidth  =  5.9058; 

% 

Stimulated  frequency  bandwidth 

sVV  60  =  VV(:,981); 
zVV_60  =  fft(sVV_60); 

% 

VV  polarized  radar  signal  at  9azimuth  =  60°  and  <peievation  =  30‘ 

sHH  60  =  HH(:,981); 
zHH_60  =  fft(sHH_60); 

% 

HH  polarized  radar  signal  at  0azimuth  =  60°  and  <peievation  =  3  0‘ 
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sHV_60  =  HV(:,981);  % 

zHV_60  =  fft(sHV_60); 


HV  polarized  radar  signal  at  Oazimuth  =  60°  and  (peievation  =30‘ 


%  Plotting  of  Radial  Range  Profile  with  VV  polarization 
figure 

plot((-length(zVV_60)/2:length(zVV_60)/2-l)*c/2/bandwidth/le9,  (abs(fftshift(zVV_60)))) 
xlabel('Range  (m)') 

ylabel('Magnitude  of  Impulse  Response') 

title('Radial  Range  Profile  for  \theta_a_z  =  60  deg  and  \theta_e_l  =  30  deg  with  VV  Polarization') 
grid  on 


%  Plotting  of  Radial  Range  Profile  with  HH  polarization 
figure 

plot((-length(zHH_60)/2:length(zHH_60)/2-l)*c/2/bandwidth/le9,  fftshift(abs((zHH_60)))) 
xlabel('Range  (m)') 

ylabel('Magnitude  of  Impulse  Response') 

title('Radial  Range  Profile  for  \theta_a_z  =  60  deg  and  \theta_e_l  =  30  deg  with  HH  Polarization') 
grid  on 


%  Plotting  of  Radial  Range  Profile  with  HV  polarization 
figure 

plot((-length(zHV_60)/2:length(zHV_60)/2-l)*c/2/bandwidth/le9,  (abs(fftshift(zHV_60)))) 
xlabel('Range  (m)') 

ylabel('Magnitude  of  Impulse  Response') 

title('Radial  Range  Profile  for  \theta_a_z  =  60  deg  and  \theta_e_l  =  30  deg  with  HV  Polarization') 
grid  on 


%  Computation  of  third  order  cumulant  and  bispectrum 
[cxVV_60]  =  cum3  ( s  V  V_6  0, (length)  s  V  V_6  0 )  - 1 ) ) ; 
bispVV_60  =  fftshift(fft2(ifftshift(cxVV_60))); 

[cxHH_60]  =  cum3  ( sHH_6  0, (length)  sHH_6  0 )  - 1 ) ) ; 
bispHH_60  =  fftshift(fft2(ifftshift(cxHH_60))); 

[cxHV_60]  =  cum3  ( sH  V_6  0, (length)  sH  V_6  0 )  - 1 ) ) ; 
bispHV_60  =  fftshift(fft2(ifftshift(cxHV_60))); 

%  Plotting  of  the  3-D  bi-range  profile  with  VV  polarization 
figure 

xVV_60  =  linspace(-length(bispVV_60)/4,length(bispVV_60)/4,length(bispVV_60))*c/2/bandwidth/le9; 
yVV_60  =  linspace(-length(bispVV_60)/4,length(bispVV_60)/4,length(bispVV_60))*c/2/bandwidth/le9; 
refVV_60  =  find(xW_60==0); 
refendVV_60  =  find(xVV_60>=2  &  xVV_60<=2.01); 

mesh(xVV_60(refVV_60:refendVV_60),yVV_60(refVV_60:refendVV_60),abs(bispVV_60(refVV_60:ref 

endVV_60,(refVV_60:refendVV_60)))); 

xlabel('Bi-range,  r_l  (m)') 

ylabel('Bi-range,  r_2  (m)') 

zlabel('Magnitude  of  Bi-range  Response') 

title('3-D  Bi-range  Profile  for  \theta_a_z  =  60  deg  and  \theta_e_l  =  30  deg  with  VV  Polarization') 

%  Plotting  of  the  3-D  bi-range  profile  with  HH  polarization 
figure 

xHH_60  =  linspace(-length(bispHH_60)/4,length(bispHH_60)/4,length(bispHH_60))*c/2/bandwidth/le9; 
yHH_60  =  linspace(-length(bispHH_60)/4,length(bispHH_60)/4,length(bispHH_60))*c/2/bandwidth/le9; 
refHH_60  =  find(xHH_60==0); 
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refendHH_60  =  find(xHH_60>=2  &  xHH_60<=2.01); 

mesh(xHH_60(refHH_60:refendHH_60),yHH_60(refHH_60:refendHH_60),abs(bispHH_60(refHH_60:ref 

endHH_60,(refHH_60:refendHH_60)))); 

xlabel('Bi-range,  r_l  (m)') 

ylabel('Bi-range,  r_2  (m)') 

zlabel('Magnitude  of  Bi-range  Response') 

title('3-D  Bi-range  Profile  for  \theta_a_z  =  60  deg  and  \theta_e_l  =  30  deg  with  HH  Polarization') 

%  Plotting  of  the  3-D  bi-range  profile  with  HV  polarization 
figure 

xHV_60  =  linspace(-length(bispHV_60)/4,length(bispHV_60)/4,length(bispHV_60))*c/2/bandwidth/le9; 
yHV_60  =  linspace(-length(bispHV_60)/4,length(bispHV_60)/4,length(bispHV_60))*c/2/bandwidth/le9; 
refHV_60  =  find(xHV_60>=-2.01  &  xHV_60<=-2.0); 
refendHV_60  =  find(xHV_60>=2.0  &  xHV_60<=2.01); 

mesh(xHV_60(refHV_60:refendHV_60),yHV_60(refHV_60:refendHV_60),abs(bispHV_60(refHV_60:ref 

endHV_60,(refHV_60:refendHV_60)))); 

xlabel('Bi-range,  r_l  (m)') 

ylabel('Bi-range,  r_2  (m)') 

zlabel('Magnitude  of  Bi-range  Response') 

title('3-D  Bi-range  Profile  for  \theta_a_z  =  60  deg  and  \theta_e_l  =  30  deg  with  HV  Polarization') 

%  Plotting  of  the  2-D  characteristic  bi-range  profile  with  VV  polarization 
figure 

contour(xVV_60(refVV_60:refendVV_60),yVV_60(refVV_60:refendVV_60),abs(bispVV_60(refVV_60:r 
efendVV_60,(refVV_60:refendVV_60))),  [i.8e5  4e5  6e5  8e5  le6  1.2e6  1.4e6  1.6e6  1.8e6]); 
xlabel('Bi-range,  r_l  (m)') 
ylabel('Bi-range,  r_2  (m)') 

title('2-D  Characteristic  Bi-range  Profile  for  \theta_a_z  =  60  deg  and  \theta_e_l  =  30  deg  with  VV 

Polarization') 

grid  on 

colorbar 


%  Plotting  of  the  2-D  characteristic  bi-range  profile  with  HH  polarization 
figure 

contour(xHH_60(refHH_60:refendHH_60),yHH_60(refHH_60:refendHH_60),abs(bispHH_60(refHH_60:r 
e fendHH_6 0 ,  (refHH_6  0 : refendHH_6 0 ) ) ) ,  [2.5e5  4e5  6e5  8e5  le6  1.2e6  1.4e6  1.6e6  1.8e6  2.0e6  2.5e6]); 
xlabel('Bi-range,  r_l  (m)') 
ylabel('Bi-range,  r_2  (m)') 

title('2-D  Characteristic  Bi-range  Profile  for  \theta_a_z  =  60  deg  and  \theta_e_J  =  30  deg  with  HH 

Polarization') 

grid  on 

colorbar 


%  Plotting  of  the  2-D  characteristic  bi-range  profile  with  HV  polarization 
figure 

contour(xHV_60(refHV_60:refendHV_60),yHV_60(refHV_60:refendHV_60),abs(bispHV_60(refHV_60:r 
efendHV_60,(refHV_60:refendHV_60))),  [3.5e6  4e6  6e6  8e6  le7  2e7  2.5e7  3e7  3.5e7]); 
xlabel('Bi-range,  r_l  (m)') 
ylabel('Bi-range,  r_2  (m)') 

title('2-D  Characteristic  Bi-range  Profile  for  \theta_a_z  =  60  deg  and  \theta_e_l  =  30  deg  with  HV 

Polarization') 

grid  on 

colorbar 
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